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ABSTRACT 


Fibrous  layers  of  the  lamina  propria  influence  dynamic  behavior  of  the 
tympanic  membrane  by  producing  strong  anisotropic  extensional  stiffness,  while 
the  mucous  and  epidermal  layers  are  primarily  responsible  for  curvature 
dependent  structural  damping  and  relatively  weak  isotropic  bending  stiffness. 
These  mechanical  properties  are  combined  with  curvilinear  shell  equilibrium 
equations  to  formulate  a  comprehensive  dynamic  continuum  model  of  the  tympanic 
membrane.  The  resulting  model  contains  several  small  parameters  that  are 
exploited  to  construct  closed  form  asymptotic  solutions  for  the  general 
problem.  Using  the  geometry  and  ultrastructure  of  the  cat  eardrum  coupled 
with  the  outer  ear,  ossicular  chain  and  tympanic  cavities,  the  asymptotic 
solution  reproduces  the  manifold  of  experimentally  observed  frequency  and 
excitation  dependent  vibrational  shapes.  Demonstration  of  the  model's  ability 
to  duplicate  experimental  results  concludes  part  I.  However,  since  the  model 
is  based  on  the  actual  geometry  and  ultrastructure,  future  investigations  may 
use  the  model  to  address  questions  related  to  tympanoplasty,  tympanosclerosis, 
tensor  tympani  cogency,  as  well  as  energy  transmission  and  tympanic  cavity 
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CHAPTER  1 


INTRODUCTION 

1.1  The  Spirit  of  Mathematical  Physics 

"Since  the  seventeenth  century,  physical  intuition  has  served  as  a 
vital  source  for  mathematical  problems  and  methods.  Recent  trends  and 
fashions  have,  however,  weakened  the  connection  between  mathematics  and 
physics;  mathemathcians,  turning  away  from  the  roots  of  mathematics  in 
intuition,  have  concentrated  on  refinement  and  emphasized  the  postulational 
side  of  mathematics,  and  at  times  have  overlooked  the  unity  of  their  science 
with  physics  and  other  fields.  In  many  cases,  physicists  have  ceased  to 
appreciate  the  attitudes  of  mathematicians.  This  rift  is  unquestionably  a 
serious  threat  to  science  as  a  whole;  the  broad  stream  of  scientific 
development  may  split  into  smaller  and  smaller  rivulets  and  dry  out.  It 
seems,  therefore,  important  to  direct  our  efforts  toward  reuniting  divergent 
trends  ..." 

R.  Courant,  1953 

1.2  Mathematical  Modeling  of  the  Eardrum 

In  the  spirit  of  mathematical  physics  described  by  R.  Courant,  a 
mathematical  model  must  have  an  intimate  connection  to  the  physical  reality 
it  is  intended  to  describe.  Dynamic  behavior  of  the  eardrum  is  governed  by 
the  conservation  laws  (mass,  momentum,  energy),  and  the  thermodynamic  entropy 
inequality.  When  applied  to  the  tympanic  membrane  (TM),  these  physical  laws 
can  be  used  to  derive  mathematical  expressions  describing  its  behavior. 
Unfortunately,  including  all  of  the  detailed  anatomy,  geometry,  and  ultra¬ 
structure  of  the  eardrum  when  deriving  the  model  is  an  immense,  and  perhaps 
imprudent,  task.  To  create  a  moderately  tractable  mathematical  model  of  the 
tympanic  membrane  it  appears  necessary  to  consider  a  simplified,  idealized, 


system.  The  degree  to  which  the  mathematics  associated  with  the  simplified 
system  describes  behavior  of  the  eardrum,  depends  upon  how  well  the  ideal¬ 
ization  represents  the  actual  physics.  The  simplification  process  and  its 
associated  reduction  in  physical  content,  raises  the  question:  What  physical 
axioms  should  be  included  in  a  tympanic  membrane  model?  Part  of  the  answer 
lies  in  how  the  model  is  to  be  used. 

Many  potential  applications  of  the  TM  model  can  be  found  in  clinical 
areas.  Most  current  clinical  methods  used  in  assesment  of  hearing  function 
rely  on  external  acoustic  stimuli  such  that  the  test  results  contain  tympanic 
membrane  factors.  A  mathematical  model  can  help  clinicians  evaluate  the 
eardrum,  and  can  also  be  used  to  remove  tympanic  membrane  factors  from  the 
data  permitting  a  more  direct  examination  of  the  inner  ear.  In  addition 
to  the  audiological  applications,  it  is  desireable  to  have  a  model  that  can 
assess  tympanosclerosis  conditions  and  impairments  caused  by  tympanoplasty. 
From  a  function  and  failure  point  of  view,  the  model  should  also  be  capable  of 
estimating  dynamic  stress  and  deflection  fields  under  a  wide  variety  of 
stimuli.  To  address  energy  transmission  and  noise  related  damage,  coupling  of 
the  TM  to  the  outer  ear,  ossicular  chain,  and  tympanic  cavity  must  be  in¬ 
cluded.  Deriving  a  model  with  these  features  requires  detailed  study  of  the 
geometry  and  ultrastructure  of  the  eardrum.  The  ability  to  describe  a  range 
of  geometries  and  fine  structures  will  give  the  model  the  flexibility  to 
describe  individual  diversity  and  differences  between  species. 

Prior  to  actual  formulation  of  the  model.  Chapter  2  reviews  histori¬ 
cally  popular  tympanic  membrane  models.  This  review  is  followed  by  two 
chapters  discussing  importance  of  various  physical  aspects  of  the  drum. 

Chapter  5  contains  the  model  derivation,  while  Chapter  6  discusses  perturba¬ 
tions  on  the  model  caused  by  the  Eustachian  tube  and  middle  ear  air  chambers. 


To  provide  completeness  and  qualitatively  describe  the  function  of  the  TM,  the 
following  section  briefly  reviews  the  macroscopic  anatomical  construction  of 
the  hearing  system. 

1.3  Macro-Anatomy  of  the  Hearing  System 

The  hearing  system  in  man  and  most  mammals  is  generally  discussed  in 
terms  of  five  broad  subregions;  the  external  ear,  the  middle  ear,  the  internal 
ear,  the  ascending  neuronal  chain,  and  the  auditory  cortex.  These  subregions 
interact  with  each  other  in  a  nearly  sequential  way.  External  acoustic  waves 
are  first  encountered  by  the  pinna  and  external  auditory  meatus  of  the  outer 
ear  as  indicated  in  Figure  1.1.  The  acoustic  signal  felt  by  the  tympanic 
membrane  differs  both  spatially  and  in  amplitude  spectrum  from  the  external 
stimuli  due  to  action  of  the  outer  ear  (Rodgers,  1981;  Stinson  and  Shaw,  1982; 


Malleus 


Fig.  1.1  Schematic  of  the  Human  Ear 


Hudde,  1982).  The  eardrum  is  generally  considered  part  of  the  middle  ear  with 
its  external  surface  bounding  on  the  external  ear  canal.  The  modified  acoustic 
signal  acting  on  the  tympanic  membrane  causes  stress  and  deflection  of  the 
membrane,  coupling  the  acoustic  signal  to  mechanical  motion.  The  stresses  in 
the  membrane  interact  directly  with  the  malleus  bone  of  the  middle  ear. 

Stresses  in  the  eardrum  at  its  connection  to  the  malleus  provide  the  force 
imparted  on  the  ossicles  by  the  eardrum. 


Mallaua  Suparior  Ligaaant 


Fig.  1.2  Schematic  of  the  Human  Middle  Ear 
In  addition  to  its  coupling  to  the  ear  canal  and  malleus,  the  tympanic 


membrane  interacts  with  the  air  (and  fluid)  within  the  primary  middle  ear  air 
cavity.  The  spatial  distribution  of  pressure  within  the  primary  cavity  of  the 
middle  ear  is  dependent  upon  the  spectrum  of  the  stimuli,  function  of  the 


Eustachian  tube,  tensor  timpani  muscle  action,  and  fluid  in  the  cavity.  For 
normal  ears  subject  to  acoustic  stimuli  well  below  5  KHZ,  the  primary  air 
chamber  acts  essentially  as  a  closed  cavity  with  homogeneous  pressure  and  den¬ 
sity.  As  shown  in  Chapter  6,  for  high  frequencies  the  pressure  distribution 
within  the  middle  ear  air  chambers  is  not  uniform. 

In  humans  the  eardrum  is  shaped  like  a  shallow  convex  cone.  The  "cone" 
is  inclined  to  the  ear  canal  cross  section  and  attached  to  the  wall  of  the 
canal  along  the  annular  ring.  Near  the  umbo  (apex  of  the  cone)  the  malleus  is 
integrally  attached  to  the  tympanic  membrane  by  woven  fibers  and  tissue. 
Attachment  of  the  malleus  to  the  tympanic  membrane  becomes  looser  towards  the 
annular  ring  (Graham,  1978). 

These  features  of  the  tympanic  membrane  allow  it  to  perform  its  primary 
function  of  converting  acoustic  stimuli  into  mechanical  force  and  motion  of 
the  ossicles.  The  first  ossicle,  the  malleus,  is  supported  by  three 
ligaments,  a  muscle,  and  the  eardrum  as  depicted  in  Figure  1.2.  The  superior, 
anterior  and  lateral  ligaments  are  passive  supports  having  relatively  low 
bending  stiffnesses  and  high  axial  stiffnesses.  In  contrast,  the  tensor  tym- 
r  i  when  stimulated  is  capable  of  forcing  the  malleus  inducing  deformation 
and  prestress  in  the  tympanic  membrane.  Deflection  of  the  malleus  relative 
to  its  size  is  very  small  such  that  the  malleus  instant  center  of  velocity 
is  essentially  stationary  at  a  given  frequency. 

Attached  almost  rigidly  to  the  malleus  is  the  incus,  the  second 
ossicle.  The  incus  is  supported  by  the  passive  posterior  ligament  in  addi¬ 
tion  to  its  strong  connection  to  the  malleus.  Despite  their  apparent  rigid 
attachment,  the  malleus-incus  joint  can  partially  decouple  under  specific 
strong  overstimulation  providing  an  apparent  internal  safety  mechanism. 


The  tip  of  the  incus  forces  the  final  ossicle,  the  stapes.  The  stapes 
is  supported  by  a  synovial  joint  at  the  incus  tip  and  is  bonded  to  the  oval 
window  membrane  of  the  cochlea.  A  second  primary  muscle  of  the  middle  ear, 
the  stapedius  muscle  forces  the  stapes  transverse  to  the  line  of  incus  action 
providing  active  control  of  input  from  the  stapes  to  the  cochlea.  When  the 
stapedius  muscle,  is  relaxed  motion  of  the  stapes  is  nearly  rectilinear, 
imparting  a  piston  like  motion  to  the  round  window.  Excitation  of  the 
stapedius  kinetically  constrains  the  stapes  to  a  rocking  motion,  reducing  the 
spatially  averaged  displacement  of  the  oval  window. 

The  cochlea  is  the  final  mechanical  element  in  the  auditory  system 
and  performs  the  conversion  from  mechanical  motion  to  a  discrete  set  of 
frequency  and  amplitude  coded  neural  signals.  The  interpretation  of  neural 
signals  generated  in  the  cochlea  is  performed  by  the  auditory  cortex.  To 
close  the  loop,  feedback  is  sent  from  the  neural  system  to  active  devices 
including  the  tensor  tympani  and  stapedius  muscles  of  the  middle  ear. 

This  brief  description  of  the  macro-anatomy  of  the  hearing  system  is 
presented  as  an  orientation  for  studying  details  of  the  mechanics  of  the 
tympanic  membrane.  Comments  and  descriptions  of  the  system  have  been  kept 
brief.  In  addition  to  the  references  mentioned,  F.  Netter  (1900)  has  pro¬ 
duced  several  very  clear  illustrations  of  the  human  hearing  system. 


CHAPTER  2 


A  CRITICAL  HISTORICAL  REVIEW  OF  TYMPANIC  MEMBRANE  MODELS 

2.1  Introduction 

Significant  efforts  to  quantitatively  describe  the  behavior  of  the 
middle  ear  in  animals  and  humans  dates  back  to  the  early  eighteen  hundreds. 
Prior  to  that  time  the  gross  anatomy  and  structure  was  known,  however  the 
description  of  the  mechanics  of  the  middle  ear  function  in  mathematical  terms 
had  not  yet  begun.  Early  theoretical  modeling  of  the  entire  middle  ear 
system,  including  the  tympanic  membrane,  was  started  essentially  without  the 
aid  of  quantitative  experimental  data  and  with  very  little  qualitative  data. 
Despite  this  lack  of  experimental  results  and  detailed  anatomical 
descriptions,  pioneer  researchers  quickly  came  to  the  (perhaps  obvious) 
conclusion  that  the  tympanic  membrane  plays  a  major  role  in  the  function  of 
the  middle  ear.  In  fact,  the  tympanic  membrane  is  responsible  for  a  large 
percentage  of  the  static  force  ratio  apparent  in  the  entire  middle  ear 
system. 

The  geometrical  shape  and  structure  of  the  eardrum  is  complex  and  has 
forced  researchers  to  limit  the  degree  of  the  actual  physics  contained  in 
mathematical  models  of  the  system.  Models  of  TM  are  created  by  considering 
a  finite  number  of  restoring,  dissipation,  and  inertia  mechanisms  coupled 
with  kinematics  of  the  deformation.  Various  models  currently  in  existence 
address  different  degrees,  or  levels,  of  physical  content.  The  content  of 
the  model  dictates  its  ability  to  reproduce  experimental  results,  and  more 
importantly,  to  predict  untested  behavior  of  the  real  physical  system. 

Models  appearing  in  the  literature  over  the  last  century  can  be  grouped  into 
five  basic  catagories  depending  on  the  type  of  physical  mechanisms  contained 
i n  the  model . 


The  groups  are: 
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1.  Lumped  Parameter  Models 

2.  Curved  Lever  Models 

3.  Stiff  Plate  Models 

4.  Membrane  Models 

5.  Shell  Bending  Models 


Even  the  simplest  lumped  parameter  models  are  able  to  describe  some  of  the 
behavior  of  the  tympanic  membrane,  however  distributed  parameter  models  are 
necessary  to  describe  actual  stresses,  vibrational  shapes,  and  effects  of 
structural  changes  on  the  system.  To  gain  a  global  understanding  of  the 
catagories,  each  type  of  model  will  briefly  be  discussed  pointing  out  some 
of  their  fundamental  strengths  and  weaknesses. 


2.2  Lumped  Parameter  Models 

Lumped  parameter  models  of  the  middle  ear  and  the  tympanic  membrane 
are  the  earliest  and  simplest  mathematical  representations  of  the  hearing 
mechanism.  Perhaps  due  to  their  simplicity,  the  lumped  (or  discrete  system) 
models  have  aged  with  the  least  controversy  and  have  been  the  outstanding 
choice  when  incorporating  eardrum  models  in  comprehensive  hearing  models  to 
date.  The  recent  emergence  of  continuous  (distributed  parameter)  systems  of 
equations  as  the  clear  leader  in  the  understanding  of  cochlear  mechanics  has 
initiated  a  move  away  from  the  discrete  models  toward  a  more  realistic  des¬ 
cription  of  the  physics  (Steele  1976;  Holmes  and  Cole  1983).  It  seems  clear 
that  detailed  descriptions  of  the  tympanic  membrane  will  also  require 
continuum  modeling  in  order  to  overcome  some  of  the  shortcomings  of  discrete 
models. 

Historically,  the  first  lumped  parameter  models  considered  the  tympanic 
membrane  to  act  like  a  piston  in  a  cylinder  forcing  a  pushrod  connected  to  the 
manubrium.  The  result  of  such  a  simple  model  is  an  algebraic  equation  relat¬ 
ing  the  average  sound  pressure  level  at  the  tympanic  membrane  to  the  force 
developed  at  the  manubrium  tip  (Esser,  1947).  More  recent  models  have 
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incorporated  average  transverse  inertia  terms  as  well  as  dissipation  terms  in 
order  to  more  accurately  reproduce  the  experimental  data.  In  many  of  the 
common  lumped  parameter  models  the  tympanic  membrane  is  represented  as  a 
second  order  ordinary  differential  equation  (or  a  system  of  n  first  order 
differential  equations).  One  of  the  more  widely  used  models  is  an  electrical 
anology  presented  by  Zwislocki  (1962)  as  outlined  in  Figure  2.1. 


I  I  I 


Block  diagram  of  the  middle-ear  mechanism. 


Figure  2.1  Zwislocki's  (1962)  circuit  model  of  the  middle  ear. 

The  mooel  has  the  ability  to  reproduce  some  experimental  data  including  com¬ 
ponents  of  reactive  and  resistive  acoustic  impedance  over  a  range  of  forcing 
frequencies.  Improvements  in  the  discrete  models  have  been  presented  by  many 
researchers  including  Zwislocki,  Moller,  Nuttall,  Fischler,  Rubenstein,  and 


Figure  2.2  Nuttall's  (1974)  modification  of  the  lumped  parameter 

ci rcui t  model . 


Stuhlman.  These  improvements  essentially  amount  to  increasing  the  degrees 
of  freedom  of  the  ordinary  differential  system.  As  an  example.  Figure  2.2 
indicates  the  lumped  parameter  system  suggested  by  A.  Nuttall  (1974)  modeling 
tympanic  muscle  effects  on  middle-ear  transfer  characteristics.  To  improve 
the  ability  of  the  lumped  parameter  models  to  reproduce  experimental  data 
above  4  KHz,  Shaw  (1977,  1982)  modeled  the  TM  as  two  coupled  pistons  doubling 
the  degrees  of  freedom  over  original  models.  The  lumped  parameter  models  are 
able  to  reproduce  some  of  the  global  features  of  the  so-called  middle  ear 
transfer  function,  however  they  have  no  hope  of  predicting  transmission 
alterations  caused  by  tympanoplasty  applied  directly  to  the  membrane  structure 
and  are  not  suited  to  describe  vibrational  modes  or  shapes  of  the  eardrum. 

The  numerical  values  of  parameters  (resistance,  inductance,  capacitance)  used 
in  the  lumped  parameter  models  are  usually  selected  on  the  basis  of  fitting 
experimental  test  data  curves.  This  method  of  selecting  the  parameter  values 
can  optimize  the  ability  of  the  model  to  reproduce  experimental  data,  but  may 
r  remove  any  real  physical  meaning  attached  to  the  numerical  values.  Without  a 

physically  based  derivation  of  the  parameter  values,  it  is  very  difficult  (or 
impossible)  to  predict  alteration  of  the  values  related  to  changes  in  the 
system.  When  using  the  discrete  models  it  is  usually  necessary  to  conduct  a 
series  of  tests  on  systems  that  have  been  altered  and  attempt  to  correlate 
changes  in  the  parameters  necessary  to  reproduce  the  results. 

In  an  effort  to  describe  the  vibrational  shape(s)  of  the  tympanic  mem¬ 
brane  in  terms  of  its  physical  features  and  sound  pressure(s)  applied,  the 
discrete  models  must  be  expanded  to  higher  dimensions.  The  obvious  step  is  to 
consider  a  series  of  continuum  models,  each  model  in  the  series  containing 
varying  degrees  of  the  actual  physics  of  eardrum.  The  earliest  of  the 
continuum  models  was  presented  by  Helmholtz,  published  in  the  middle  eighteen 


hundreds. 
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2.3  Curved  Lever  Models 

The  pioneer  continuum  type  model  of  the  tympanic  membrane  was  presented 
by  Helmholtz  in  a  paper  entitled  "The  Mechanism  of  the  Ossicles  of  the  Ear  and 
Membrana  Tympani"  in  which  he  presented  the  'curved  membrane'  theory  of  the 
mechanism  of  the  tympanic  membrane.  No  explicit  model  was  presented,  in  fact 
a  complete  mathematical  representation  and  solution  of  the  mechanism  remains 
open;  however,  a  physical  mechanism  describing  action  of  the  ear  drum  was 
given.  The  theory  rests  on  the  convex  cone  shape  of  the  tympanic  membrane. 

If  we  imagine  that  the  membrane  consists  of  a  number  of  curved  radial  fibers 
connecting  the  malleus  to  the  annular  ring,  then  it  is  possible  to  describe  a 
'lever  type'  action  of  each  of  the  radial  fibers  acting  on  the  malleus.  To 
illustrate  this.  Figure  2.3  shows  a  single  radial  curved  fiber  attached 
rigidly  to  the  annular  ring  exposed  to  a  pressure  differential  P.  If  the 
fiber  is  essentially  inextensible  then  inward  displacements  of  the  manubrium 
will  be  accompanied  with  an  increased  radius  of  curvature  of  the  fiber,  and 
outward  displacements  of  the  manubrium  will  be  accompanied  with  a  diminished 
radius  of  curvature.  These  changes  in  curvature  alter  the  distance  between 
the  endpoints  of  particular  radial  fibers  causing  the  manubrium  to  move 
accordingly  (Fig.  2.3).  The  kinematics  described  converts  large  displacement 
(low  force)  motions  of  the  membrane  into  small  displacement  (high  force) 
motions  of  the  manubrium  and  hence  inherits  the  name  curved  lever  model. 

The  original  curved  membrane  theory  presented  by  Helmholtz  assumed  that 
the  radial  fibers  were  inextensible  and  that  changes  in  curvature  of  the 
radial  fibers  was  opposed  by  elastic  behavior  of  the  circumferential  fibers. 
This  premise  of  Helmholtz's  curved  membrane  theory  is  reasonable  only  if  the 
radial  stiffness  is  much  greater  than  the  circumferential  stiffness. 
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Experimental  data  indicates  that  the  bending  modulus  of  the  drum  material  in 
the  radial  direction  is  the  same  order  of  magnitude  as  the  bending  modulus  in 
the  circumferential  direction  (von  Bekesy,  1948;  Kirikae,  1960).  Due  to  the 
fiber  composite  structure,  this  bending  stiffness  observation  does  not  imply 
the  same  conclusion  for  membrane  type  stiffnesses  (see  Chapter  3).  The 
membrane  stiffnesses  may  be  more  orthotropic  than  the  bending  stiffnesses. 
Several  authors  have  argued  that  the  relatively  isotropic  bending  stiffness 
discounts  the  curved  membrane  theory  (Funnel  1  and  Laszlo,  1978;  Khanna  and 
Tonndorf,  1972).  It  should  be  pointed  out  that  axial  stiffnesses  are 
appropriate  to  use  in  such  an  arguement,  not  bending  stiffnesses.  Secondly, 
the  restoring  mechanism  in  the  curved  lever  theory  need  not  be  restricted  to 
extension  of  the  circular  fibers,  but  is  more  likely  connected  to  the  total 
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Figure  2.3.  Using  a  static  force  balance,  and  assuming  small  deformations, 
the  force  per  unit  length  on  the  mallus  due  to  one  radial 
fiber  has  the  form  F  ”  P[cos( a)  ( l  .u.°.)~  +  z/2]. 


deformation  through  a  general  anisotropic  constitutive  law.  Hence,  histor¬ 
ical  arguments  used  to  justify  neglecting  radial  curvature  and  anisotropic 
construction  are  not  valid  these  items  may  indeed  play  a  first  order  role. 


In  the  middle  1940's,  prior  to  Urn's  work,  Esser  (1947)  formalized 
Helmholtz's  curved  membrane  theory  into  an  axi symmetric  set  of  equilibrium 
equations  relating  the  stress  to  the  curvature,  thickness,  and  applied  pres¬ 
sure.  The  static  equations  used  by  Esser  are  axi symmetric  pressure  vessel 
equations  and  are  given  by: 


—  (rhoj  -  ho*  *  0 

Hr  IN  t 


(2.3.1) 


ho  d  (cos  a)  +  ho*  1  cos  a  =  P 
m  clr  t  7 

where  P  is  the  applied  pressure,  at  is  the  tangential  stress,  om  is  the 
meridian  stress,  h  is  the  thickness  a  is  the  local  angle  of  a  normal  from  the 
membrane  surface,  and  r  is  the  radial  position  coordinate. 

Since  the  stresses  depend  on  the  deformation,  and  since  a  also  depends 
on  the  deformation,  the  equations  are  nonlinear.  Esser  considered  several 
simplified  forms  as  assumed  stress  states  in  order  to  estimate  the  effective 
static  "lever  ratio"  of  the  membrane  as  predicted  by  Helmholtz.  Essers' 
approximate  solutions  indicate  that  curvature  may  indeed  play  a  role  in  the 
function  of  the  tympanic  membrane.  Although  the  results  of  Esser  may  be 
written  in  terms  of  a  "lever  ratio",  it  should  be  noted  that  dynamic  effects, 
quasilinearities  in  stress,  nonsymmetric  shape,  and  kinematic  coupling  were 
not  considered.  (For  an  indication  of  the  mathematics  involved  in  expanding 
the  equations  to  the  case  of  "Thin  Elastic  Membranes"  see  Pujara  and  Ladner, 
1978,  or  Yang  and  Feng,  1970.)  Omission  of  these  factors  makes  the  "model" 
inadequate  for  comprehensive  study,  and  should  simply  be  viewed  as  evidence 
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of  the  importance  of  curvature  and  membrane  type  stress. 

The  curved  membrane  theory  was  rekindled  by  Tonndorf  and  Khanna  (1972) 
with  their  suggestion  that  the  restoring  force  may  be  due  to  shearing  of  the 
radial  fibers  over  the  circumferential  and  parabolic  fibers  rather  than  the 
original  extensional  restoring  force.  The  modified  theory  maintains  the 
"lever  type"  action  of  the  membrane  and  allows  both  sets  of  primary  fibers  to 
have  similar  elastic  moduli  imbedded  in  a  base  material  having  a  much  lower 
stiffness.  One  difficulty  with  the  modified  theory  is  that  the  circular 
fibers  must  slide  a  distance  of  approximately  two  diameters  over  the  radial 
fibers  during  a  typical  deformation  if  both  sets  of  fibers  are  indeed  inex- 
tensible  (Tonndorf  and  Khanna,  1972).  This  distance  seems  large  when  one 
considers  that  almost  no  material  exists  between  the  two  sets  of  fibers  (Lim, 
1968  a,b,  1970).  Calculations  done  in  Chapter  3  indicate  that  the  membrane 
indeed  contains  an  internal  lever  system,  however,  the  restoring  mechanism  may 
be  dominated  by  tension  and  bending  terms  rather  than  the  shear  or  extensional 
mechanisms  mentioned  above.  The  experimentally  based  results  of  Chapter  3 
also  show  that  the  membrane  stress  induced  in  the  radial  fibers  by  a  positive 
pressure  is  extensional  rather  than  compressive,  discounting  the  original 
curved  lever  mechanism  postulated  by  Helmholtz. 

2.4  Membrane  Models 

The  "curved  membrane"  model  initialized  by  Helmholtz  includes  kinema¬ 
tics  associated  with  curvature  and  normal  membrane  type  stresses.  Elementary 
stresses  in  the  plane  of  a  thin  material  surface  are  considered  of  the  mem¬ 
brane  type  if  their  resultant  bending  moment  is  zero  on  every  normal  cut 


through  the  surface;  mathematically  this  means: 

h/2 

/  /  oz  dz  ds  =  0  (2.4.1) 

S  -h/2 

where  o  is  the  stress  tangent  to  the  middle  surface,  z  is  the  distance  from 
the  neutral  surface,  h  is  the  thickness,  and  s  is  an  arbitrary  curve  in  the 
middle  surface.  With  this  definition,  a  mathematical  model  is  considered  a 
membrane  model  if  membrane  type  stresses  are  the  only  stresses  considered  when 
formulating  the  model.  Simple  membrane  models  assume  that  the  membrane  stress 
is  uniform  through  the  thickness,  and  that  variations  in  the  stress  resulting 
from  deformation  are  small  relative  to  some  initial  stress.  More  detailed 
membrane  theories,  such  as  presented  by  Yang  and  Feng  (1970),  are  able  to 
describe  changes  in  membrane  stress  due  to  applied  pressures  and  resulting 
deformation.  Membrane  equations  or  theories  used  to  model  the  tympanic 
membrane  have  addressed  simple  linear  membrane  theory  for  "small"  deformations 
and  "large"  initial  prestress. 

Early  attempts  to  model  the  tympanic  membrane  as  a  vibrating  membrane, 
having  only  membrane  type  stress  as  the  restoring  mechanism,  are  represented 
by  Franks'  (1923)  flat  membrane  model.  The  differential  equation  considered 
by  Frank  was  the  simple  linear  membrane  equation  subject  to  a  fixed  circular 
boundary  condition  and  radial  manubrium  constraint 

chv2u  +  P  =  hp  ,  (2.42) 

3t2 

where  u  is  deformation  normal  to  the  membrane  surface,  p  is  the  mass  density, 
o  is  the  in-plane  tension,  P  is  the  applied  normal  pressure,  h  is  the  thick¬ 
ness,  t  is  time  and  v2  is  the  Laplace  operator. 

A  numerical  solution  of  the  curved  membrane  problem,  subject  to  the 
relatively  complex  boundary  conditions  imposed  at  the  manubrium  and  annular 
ring,  was  not  practical  at  the  time,  forcing  Frank  to  simplify  the  problem 


such  that  an  analytical  solution  could  be  found.  For  simplicity  Frank  elected 
to  consider  a  round,  flat,  isotropic  membrane  having  the  manubrium  attached 
along  one  radial  segment  pivoting  at  the  annular  ring.  His  results  can  be 
used  to  suggest  areas  of  maximum  deflection  and  estimate  possible  vibrational 
shapes  at  low  forcing  frequencies;  the  restrictive  assumptions  however, 
prevent  the  model  from  being  of  much  practical  use. 

Proceeding  along  the  same  lines  as  Frank,  Funnel!  (1975)  relaxed 
assumptions  on  the  circular  shape  and  manubrium  attachment  in  order  to  model 
the  cat  eardrum  (as  well  as  several  other  mammals).  Funnel  1  also  assumed  a 
flat  isotropic  membrane  operator,  however  he  imposed  more  realistic  boundary 
conditions  at  the  malleus  and  annular  ring.  A  solution  of  the  membrane 
equations  subject  to  boundary  conditions  was  estimated  using  a  finite  element 


600  Hz  4001  Hz  5176  Hz 


Figure  2.4  Experimentally  recorded  time  averaged  vibration  shapes  of  the 
the  cat  eardrum  (retouched  from  Khanna  and  Tonndorf,  1972). 

approximation.  The  tension  used  in  the  model  was  adjusted  such  that  the  cal¬ 
culated  vibrational  shapes  corresponded  to  the  shapes  observed  in  experiment 


at  selected  frequencies.  Using  this  method.  Funnel  1  was  able  to  reproduce  a 
transition  from  the  fundamental  vibratory  shape  into  a  rather  complex  shape 
similar  to  that  observed  experimentally  (Fig.  2.4).  The  amplitude  and  trans¬ 
mitted  force  predicted  by  the  plane  membrane  model  corresponded  in  a  gross 
qualitative  way  to  experimental  results,  however,  the  connection  "was  not  of 
the  quality  desired  in  a  tympanic  membrane  model"  (Funnell,  1975).  The 
vibrational  shapes  corresponded  quite  well  to  the  observed  shapes;  it  should 
be  noted  however,  that  the  shapes  are  generally  insensitive  to  the  particular 
model  and  restoring  forces  selected.  The  qualitative  shape  of  vibration  at 
low  frequency  is  influenced  strongly  by  the  boundary  shape  and  constraints, 
and  is  less  sensitive  to  the  type  of  restoring  forces. 

The  qualitative  picture  obtained  using  plane  membrane  models  is  quite 
informative,  but  we  must  remember  that  the  model  neglects  initial  curvature  of 
the  drum,  ignores  bending  stiffness,  and  disregards  bending  dissipation.  In 
addition  to  these  primary  factors,  the  membrane  type  models  discussed  neglect 
quasilinear  stress  terms,  anisotropic  construction,  rotary  inertia,  and  simi¬ 
lar  terms  that  could  be  included  in  the  membrane  type  equation(s). 

2.5  Shell  Bending  Models 

In  addition  to  membrane  type  restoring  forces,  it  is  clear  from  exper¬ 
imental  observation  that  bending  terms  also  play  a  role  in  deformation  of  the 
tympanic  membrane.  The  fact  that  the  eardrum  maintains  its  shape  after  it 
has  been  removed  from  its  supporting  structure  can  not  be  explained  using 
membrane  stresses  alone,  and  requires  inclusion  of  bending  type  restoring 
forces  (Gram,  1968).  The  bending  mechanism,  however,  may  be  secondary  in  the 
vibrating  drum  due  to  pressure  induced  membrane  stress  and  tensor  tympani  pre¬ 


stress. 
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If  a  surface  with  thickness  undergoes  a  change  in  curvature,  then  in 
general  nonuniform  strain  is  induced  across  the  thickness.  For  materials 
having  some  elasticity,  the  strain  will  be  resisted  by  stress  such  that 
strain  energy  will  be  stored  or  released  due  to  the  changes  in  curvature. 

If  a  particular  model  contains  terms  describing  this  energy,  then  the  model  is 
said  to  contain  "bending  type"  restoring  forces.  For  example,  the  Kirchhoff 
plate  equation  given  below  is  the  simpliest  of  the  plate  bending  theories 

-  D  V2V2u  +  P  =  ph  32u  ,  (2.5.1) 

3t2 

where  D  is  the  bending  stiffness,  u  is  the  deflection  normal  to  the  plate,  P 
is  the  applied  pressure,  h  is  the  thickness,  p  is  the  mass  density,  t  is  time 
and  V2  is  the  Laplace  operator.  In  the  form  given,  the  Kirchhoff  plate  equa¬ 
tion  assumes  the  plate  is  flat,  isotropic,  uniform,  undergoes  small  amplitude 
deformations,  has  no  membrane  type  stresses  or  deformations,  and  moves  at  low 
frequencies.  Some  of  these  assumptions  can  be  relaxed  by  considering  the  Von 
Karman  plate  or  shell  equations  containing  bending  and  membrane  type  restoring 
forces.  Early  models  of  the  tympanic  membrane  have  neglected  one  or  the  other 
of  the  two  primary  restoring  mechanisms.  Some  membrane  type  models  have 
already  been  discussed.  In  addition  to  the  models  containing  purely  membrane 
type  stresses,  models  have  been  formulated  containing  purely  bending  type 
terms.  In  reality  both  effects  exist  in  the  tympanic  membrane,  but  for  the 
sake  of  mathematical  simplicity,  researchers  have  historically  considered  the 
most  degenerate  forms  that  seemed  reasonable.  Following  experimental  obser¬ 
vations  by  Bekesy  that  indicated  larger  bending  effects  than  had  previously 
been  considered.  Gram  (1968)  modeled  the  eardrum  as  a  simple  flat,  isotropic 
Kirchhoff  plate.  He  computed  an  approximate  solution  subject  to  simply 
supported  edge  conditions  and  simple  support  at  the  manubrium  using  the 
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Ray! eigh-Ri tz  method.  From  the  results  of  his  model,  he  calculated  the 
stiffness  of  the  plate  that  would  be  required  in  order  for  the  model  results 
to  match  the  experimental  impedance  data.  In  order  to  match  the  data,  the 
model  plate  required  a  stiffness  two  orders  of  magnitude  larger  than  the 
actual  stiffness  measured  in  experiment  (Bekesy,  1948;  Kirikae,  1960).  This 
discrepency  strongly  suggests  that  the  model  is  not  appropriate  for  the  ear¬ 
drum. 

The  first  attempt  to  include  both  the  "curved  cone"  shape  of  the  drum 
along  with  bending  type  restoring  forces  was  presented  by  Funnel  1  and  Laszlo 
(1977)  in  the  form  of  a  finite  element  approximate  solution  of  a  thin  shell 
model  of  the  cat  eardrum.  The  finite  element  method  is  clearly  an  effective 
way  for  estimating  solutions  to  some  particular  differential  system  of  equa¬ 
tions  and  constraints,  however  the  ability  of  the  approximate  solution  to 
describe  the  behavior  of  the  physical  system  under  study  depends  on  the 
original  set  of  differential  or  integral  equations  used.  Funnell  and  Laszlo 
elected  to  use  a  static  (or  low  frequency)  thin  shell  functional  in  which 
shell  bending  stresses  were  the  primary  restoring  forces  considered  in  the 
original  equations.  The  functional  used  contained  membrane  type  terms, 
however,  Funnell  and  Laszlo  did  not  impose  any  prestress  and  did  not  consider 
the  fibrous  ultrastructure  such  that  the  membrane  restoring  forces  were  of 
second  order.  The  use  of  bending  stress  as  the  physical  mechanism  storing 
strain  energy  during  the  deformation  had  often  been  considered  secondary  to 
other  mechanisms  in  earlier  models.  Helmholtz's  original  curved  membrane 
theory  postulated  tension  in  the  circumferential  fibers  as  the  dominant 
restoring  force,  and  Tonndorf  and  Khanna's  version  of  the  curved  membrane 
theory  considers  internal  axial  shear  as  the  dominant  storage  mechanism.  Thi 
change  in  the  model  content  is  somewhat  masked  by  the  fact  that  Funnell  and 


taszlo  were  the  first  to  apply  a  computationally  intensive  method  (finite 
element)  to  the  eardrum  problem. 

The  restoring  force(s)  considered  in  a  particular  model  alter  the 
mechanism  by  which  the  tympanic  membrane  interacts  with  the  manubrium  and  the 
annular  ring,  and  thus  may  alter  predicted  transmission  characteristics  of  the 
drum.  The  curved  membrane  models  predict  the  tympanic  membrane  will  force  the 
manubrium  through  a  membrane  stress  in  the  radial  fibers  while  the  shell 
bending  model  forces  the  manubrium  through  transverse  shear  and  moment  terms. 
This  difference  in  the  mechanism  of  force  transmission  may  be  the  primary 


Figure  2.5  Vibration  pattern  calcuated  by  Funnell  and  Laszlo,  1977.  The 

contour  lines  are  lines  of  constant  vibration  amplitude,  equally 
spaced  on  an  amplitude  scale.  Only  displacements  shown  corres¬ 
pond  to  the  component  of  displacement  measured  in  the  holographic 
experiments  of  Khanna  and  Tonndorf.  The  small  triangle  repre¬ 
sents  the  point  of  maximal  vibration  amplitude,  correspondi ng  to 
0.32  urn  for  a  low-frequency  pure-tone  input  of  100  dB  SPL.  (From 
Funnell  and  Laszlo,  1977.) 


difference  between  the  model  types.  In  a  purely  small  deformation  bending 
model  it  is  not  possible  to  have  the  internal  'lever  action'  of  the  membrane 
that  is  inherently  contained  in  membrane  type  models.  Although  the  force 
transmission  values  are  expected  to  differ  considerably  when  the  restoring 
mechanism  is  changed,  the  fundamental  vibrational  mode  shape  will  qualitatively 
remain  unchanged.  Figure.  2.5  shows  a  numerical  plot  of  deformation  of  the 
tympanic  membrane  of  a  cat  when  forced  at  low  frequencies  (Funnel  1  and  Laszlo, 
1977). 

An  apparent  shortcoming  of  the  isotropic  thin  shell  model  is  that 
natural  asymmetric  modes  of  vibration  (realized  as  peaks  on  the  "Malleus  tip 
vs.  TM  amplitude  graph")  appear  with  spacing  of  (2n)2  (n  =  1,2,3).  The  (2n)2 
spacing  is  wider  than  that  seen  experimental ly  (Manley  and  Johnstone,  1974). 

The  thin  shell  model  with  no  damping  is  also  unable  to  reproduce  decreases 
in  malleus  motion  with  frequency  and  average  phase  shift  observed  experiment¬ 
ally  (Manley  and  Johnstone,  1974;  Moller,  1972). 

Despite  some  inconsistencies  with  experimental  work,  the  vibrational 
shape  calculated  by  Funnell  and  Laszlo  is  very  similar  to  the  shape  that  has 
been  observed  experimentally  for  low  forcing  frequencies.  At  higher  frequen¬ 
cies,  when  inertia  is  added  the  bending  shell  model  does  not  reproduce  the 
complicated  vibrational  shapes  observed  in  vibrating  eardrums.  It  is  expected 
that  this  behavior  is  directly  tied  to  the  anisotropic  membrane  stress  and 
construction  of  the  drum,  and  cannot  be  reproduced  using  a  simple  isotropic 
shell  bending  model . 

2.6  Stiff  Plate  Models 

In  the  early  1940's  George  von  Bekesy  recorded  the  first  measurements  of 
the  vibrating  human  eardrum  using  a  capacitive  probe  on  fresh  cadaver  specimens. 


His  measurements,  at  frequencies  below  2000  Hz,  indicated  that  the  tympanic 
membrane  has  a  vibrational  shape  similar  to  a  flat  plate  rotating  around  a 
pivot  near  the  head  of  the  malleus.  Displacement  of  the  eardrum  in  the  stiff 
plate  mode  is  sketched  in  Figure  2.6.  The  shape  von  Bekesy  inferred  from  his 
data  is  now  known  to  be  incorrect,  however  it  is  worth  mentioning  due  to  its 
historical  significance  and  influence  on  middle  ear  research. 
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Figure  2.6  Stiff  plate  vibration  pattern  of  the  human  tympanic  membrane 

as  reported  by  von  Bekesy  (1941). 

From  a  modeling  point  of  view,  the  stiff  plate  concept,  if  correct, 
would  mean  the  tympanic  membrane  would  admit  nicely  to  a  discrete  or  lumped 
parameter  model.  In  fact,  belief  that  the  drum  vibrates  as  a  stiff  plate 
gave  rise  to  simple  discrete  modeling  concepts  like  the  drum  effective  area, 
as  a  fraction  of  the  actual  frontal  area,  to  be  used  as  a  hydraulic  lever 
area.  Existence  of  the  shape  suggested  by  von  Bekesy  requires  the  eardrum  to 
bend  and  deform  drastically  in  the  lower  part  and  has  directed  past  research¬ 
ers  to  look  for  a  'lower  fold'  that  is  now  known  not  to  exist.  If  such  a 
fold  did  exist  it  would  place  significant  structural  constraints  on  the 
tympanic  membrane.  Inspection  of  the  membrane  structure  itself  would  suggest 
that  a  'lower  fold'  is  not  expected.  It  is  assumed  that  von  Bekesy  arrived 


at  the  stiff  plate  theory  on  the  basis  of  experimental  data  that  was  not  as 
point  sensitive  as  he  expected. 

2.7  Concluding  Remarks 

Of  the  four  model  types  discussed,  the  lumped  parameter  models  have 
commanded  the  vast  majority  of  attention  in  the  literature  and  have  provided 
tools  for  interpreting  some  clinical  data.  The  lumped  models,  however,  are  by 
nature  unable  to  describe  detailed  vibrational  behavior  of  the  drum  and  are 
often  linked  to  the  actual  physics  only  through  parameter  curve  fitting.  In 
the  interest  of  gaining  a  more  complete  understanding  of  the  middle  ear 
function,  higher  order  discrete  systems  or  continuous  system  models  must  be 
addressed.  Mathematically,  the  only  complete  continuum  formulation  to  date  is 
the  bending  shell  model  presented  by  Funnel!  and  Laszlo  that  uses  standard 
small  deformation  isotropic  shell  theory.  The  highly  anisotropic  fiber 
construction  of  the  drum  has  motivated  others  to  formulate  partial  models 
based  on  restoring  mechanisms  other  than  the  shell  bending  device.  Complete 
mathematical  representations  of  the  tension  of  shear  (other  than  bending) 
models  has  not  been  presented  and,  hence,  results  from  such  models  are  not 
available  for  evaluation. 

One  important  point  that  appears  to  be  lacking  in  the  literature  is  a 
systematic  evaluation  of  possible  restoring  mechanisms  along  with  inertial  and 
damping  devices  within  the  actual  tympanic  membrane.  Analysis  and  discussion 
contained  in  the  following  four  chapters  will  center  around  determining  the 
primary  mechanisms  governing  action  of  the  tympanic  membrane  and  deriving 
mathematical  descriptions  of  the  physics.  As  a  starting  point,  the  ultra- 
structure  of  the  eardrum  is  studied  in  order  to  estimate  its  mechanical  pro¬ 
perties.  Once  the  properties  are  known,  restoring  and  dissipation  is 


addressed  followed  by  formulation  of  dynamic  equilibrium  equations  to  complete 
the  model.  A  special  closed  form  asymptotic  solution  of  the  model  is  deter¬ 
mined  in  Chapter  5  by  using  natural  small  parameters  appearing  in  the  system 
equations. 


CHAPTER  3 


TYMPANIC  MEMBRANE  ULTRASTRUCTURE  AND  STIFFNESS  DECOMPOSITION 
3.1  Introduction 

The  human  tympanic  membrane  can  be  divided  into  two  regions  having 
distinct  tissue  ultrastructures.  Pars  tensa  represents  a  majority  of  the  drum 
with  the  exception  of  a  small  triangular  region  superior  to  the  short  process 
of  the  malleus  that  consists  of  pars  flaccida  (Fig.  3.1).  The  annular 
ligament  provides  a  distinct  border  between  the  external  auditory  meatus  and 
the  pars  tensa.  Bordering  the  pars  flaccida  however,  the  annular  ligament 
terminates  making  the  pars  flaccida  integral  with  the  external  auditory 
meatal  skin  (Lim,  1970).  For  adults  the  major  diameter  of  the  annular  ring 
ranges  from  9.0  to  10.2  mm  and  the  minor  diameter  ranges  from  8.5  to  9.0  mm. 


Figure  3.1  Sketch  showing  gross  anatomical  regions  of  the  human  tympanic 

membrane. 


Tympanic  membranes  of  other  mammals  including  cats,  guinea  pigs, 
rabbits,  and  squirrel  monkeys  have  gross  anatomical  properties  similar  to  that 
of  human,  although  the  size  and  shape  of  the  membranes  vary  as  well  as  the 
fraction  of  total  area  accounted  for  by  pars  flaccida  (Funnel  1  ,  1975).  The 
shape  of  the  pars  flaccida  region  is  generally  triangular  or  half  elliptical, 
depending  on  the  particular  mammal,  with  the  base  of  the  triangle  or  half 
ellipse  contacting  the  skin  of  the  outer  ear  canal.  The  pars  tensa  of  the 
rabbit  has  approximately  four  times  the  surface  area  as  the  pars  flaccida, 
approximately  eight  times  for  the  cat,  and  fifteen  times  for  the  squirrel 
monkey.  An  exception  reported  by  Lim  (1968b)  is  the  sheep  having  a  pars  tensa 
of  approximately  twice  the  area  of  the  pars  flaccida. 

Correcting  some  historical  misconceptions,  Lim  (1968a,  1968b,  1970)  has 
shown  by  light  and  electron  microscope  investigations  that  the  pars  flaccida 
is  generally  thicker  than  the  pars  tensa.  The  pars  flaccida  of  the  rabbit  is 
20  to  30  times  thicker  than  the  pars  tensa  and  the  sheep  pars  flaccida  is  10 
to  20  times  thicker  than  the  pars  tensa  (Lim,  1968b).  Historically,  the 
thickness  of  the  pars  flaccida  was  believed  thinner  based  on  experimentally 
observed  large  static  deflections  of  the  pars  flaccida  (Shrapnell,  1832). 

This  misinterpretation  of  experimental  observation  resulted  from  not  knowing 
or  understanding  mechanical  stiffness  properties  of  the  membrane  tissues.  It 
was  later  discovered  that  the  ultrastructure  of  pars  flaccida  and  pars  tensa 
regions  are  quite  different  resulting  in  different  mechanical  properties 
(Nishiyama,  1937;  Lim,  1968a ,b). 

From  a  modeling  point  of  view,  if  we  wish  to  accurately  predict  defor¬ 
mation  and  behavior  of  the  tympanic  membrane  the  mechanical  properties  must  be 
known.  Viewing  the  eardrum  as  a  general  Von  Karman  composite  shell,  there  may 
exist  as  many  as  36  independent  mechanical  stiffness  constants  at  each  point 


in  the  drum.  Due  to  the  membrane's  small  size,  fragile  cellular  structure, 
and  post-mortem  changes,  current  testing  technology  is  not  able  to  measure  all 
of  the  mechanical  stiffnesses.  Even  if  such  an  ability  did  exist,  individual 
diversity  of  ul trastructure  including  scar  tissue,  tympanoplasty,  and  normal 
differences  would  also  need  to  be  measured  and  understood.  Rather  than 
attempting  to  develop  an  advanced  mechanical  testing  ability,  analytical 
estimation  of  the  mechanical  properties  based  on  the  local  membrane  ultra¬ 
structure  seems  to  be  a  much  more  tractable  and  reasonable  task.  Study  of  the 
fine  ultrastructure  of  the  pars  tensa  and  pars  flaccida  allows  estimation  of 
their  mechanical  stiffnesses  based  on  properties  and  geometry  of  the  tissue 
components. 

To  derive  estimates  of  mechanical  properties  of  the  tympanic  membrane 
at  various  locations,  we  must  study  its  ultrastructure.  Figure  3.2  shows  a 
cross  sectional  view  of  the  pars  tensa  of  a  squirrel  monkey  (Lim  1968a). 
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Figure  3.2  A  schematic  of  a  radial  cross  sectional  view  of  the  pars  tensa 
of  the  squirrel  monkey  (from  Lim,  1968a). 


The  pars  tensa  region  of  the  tympanic  membrane  is  constructed  of  several 
layers  of  tissue  bonded  together  in  a  laminate  structure.  The  three  primary 
layers  are  the  outer  epidermal  layer,  the  middle  lamina  propria,  and  the  inner 
mucous  layer.  (The  following  descriptions  are  condensed  from  Urn,  1968a.) 
Epidermal  layer 

The  epidermis  is  divided  into  the  stratum  corneum,  the  stratum 
granulosum,  and  the  stratum  Malpighii.  The  stratum  corneum  is  a  thin  layer 
consisting  of  degenerate  compressed  cellular  structures.  Directly  inside  of 
the  stratum  corneum,  two  layers  of  keratohyalin  granules  form  the  stratum 
granulasum.  The  outer  part  of  the  stratum  Malphigii  is  stratum  spi nosum 
consisting  of  several  layers  of  squamous  cells.  Cells  in  the  stratum  spi nosum 
are  not  as  compressed  as  in  the  stratum  granulosum  having  more  intercellular 
space.  The  squamous  cells  are  bordered  by  abundant  desmosomes  which  Lim 
(1968a)  suggests  act  as  "tightening  apparatus  which  bound  each  cell  together". 
The  stratum  basal e  is  the  innermost  layer  of  the  epidermis  and  consists  of  a 
single  layer  of  cells  separated  from  the  connective  tissue  by  a  membrane. 
Lamina  propria 

The  lamina  propria  contains  the  strong  fibers  of  the  tympanic  membrane 
separated  from  the  epidermus  and  mucous  by  connective  tissues.  The  subepi- 
dermal  connective  layer  consists  of  loose  typical  collagen  fibrils,  loose 
protocollagen  fibrils,  fibroblasts,  and  nerve  fibers.  Neurosecretory 
granules,  neurofilaments,  and  neurovesicles  exist  in  the  nerve  fibers.  Tissue 
of  the  connective  layer  blends  between  radial  fibers  of  the  outer  radiate 
collageneous  layer.  Between  the  radial  and  circular  fiber  layers,  a  thin 
parabolic  fiber  layer  exists  in  humans  (Lim,  1970).  The  subepidermal 
connective  layer  and  the  radial  fiber  layer  are  generally  more  developed  than 


the  parabolic  and  circular  fiber  layers.  The  radial  and  circumferential  fibers 
consist  of  fine  fibrils  having  rectangular  or  square  cross  sections.  The  fine 
fibrils  are  mixed  with  typical  collagen  fibrils.  In  general,  the  fine  rectan¬ 
gular  fibrils  dominate  in  both  the  radial  and  circular  layers,  however,  in 
humans  more  collagen  fibrils  are  found  in  the  circular  layer.  The  rectangular 
fibrils  do  not  show  tropocollagen  banding  and  are  finer  in  size  (approx.  50A 
vs.  200  -  400A  for  typical  collagen).  Some  indirect  evidence  suggests  the 
fine  fibrils  are  reticulin  based  (Ciges,  1965)  or  keratin  based  (Hamiton,  1967) 
The  circular  fiber  layer  borders  on  the  innermost  lamina  submucosal  connective 
tissue  layer.  Content  and  structure  of  the  submucosal  layer  is  similar  to  the 
subepidermal  layer,  except  that  the  submucosal  layer  is  thinner. 

Mucous  Layer 

The  innermost  layer  bordering  on  the  tympanic  cavity  is  the  mucous 
layer  composed  mostly  of  simple  squamous  epithelial  cells.  Near  the  annular 
ring  some  cuboidal  cells  and  cilia  are  found.  No  tonofi laments,  kerato  hyalin 
granules  or  secretory  granules  are  found  in  this  thin  layer. 

These  comments  on  the  ultrastructure  of  the  pars  tensa  are  very  con¬ 
densed.  Quantitative  content  and  size  of  the  various  layers  changes  with 
position  on  the  tympanic  membrane  and  with  species.  In  addition  to  local 
structure,  blood  vessel  (Politzer,  1892)  and  nerve  fiber  (Wilson,  1907) 
networks  have  been  documented  macroscopically.  More  detailed  descriptions  of 
the  cellular  structure  including  brief  historical  accounts  and  photomicrographs 
can  be  found  in  Lim  (1968a, b,  1970). 

The  ultrastructure  of  the  pars  flaccida  differs  from  that  of  the  pars 
tensa  primarily  in  the  lamina  propria.  Lim  (1968b)  has  shown  that  a  lamina 
propria  does  exist  in  the  pars  flaccida,  however,  contains  no  organized  fiber 


structures  similar  to  the  radial  or  circular  fibers  found  in  the  pars  tensa. 
The  fine  "rectangular"  fibrils  of  the  radial  and  circular  fibers  are  not  found 
in  the  pars  flaccida.  Larger  collagenous  fibers  are  organized  into  a  thin 
cloth  within  the  lamina  propria  of  the  pars  flaccida  imbedded  in  a  substantial 
amount  of  extracellular  amorphous  ground  substance.  This  ground  substance  and 
lack  of  a  middle  fibrous  layer  make  the  thicker  pars  flaccida  more  translucent 
than  the  pars  tensa.  The  lamina  propria  of  the  pars  flaccida  also  contains 
abundant  elastic  fibers  not  found  in  the  pars  tensa. 

Both  the  pars  tensa  and  pars  flaccida  consist  of  relatively  stiff 
fibers  imbedded  in  relatively  flexible  base  materials.  The  geometry,  orienta¬ 
tion,  and  fiber  fractions  dictate  mechanical  behavior  of  the  material (s).  In 
general,  for  a  fiber  composite  material  of  this  type,  the  first  order  linear 


constitutive  law  has  the  form 


N  *  Q  e  +  B  k 


M  =  B  e  +  D  k  , 


(3.1.1) 


where  N  is  the  normal  force  vector,  M  is  the  moment  vector,  J)  is  the  membrane 
stiffness  matrix,  £  is  the  membrane-bending  coupling  matrix,  D.  is  the  bending 
stiffness  matrix,  e  is  the  strain  vector,  and  k  is  the  curvature  vector. 

We  would  like  to  determine  constants  appearing  in  the  matricies  in 
terms  of  the  ul trastructure  of  the  tympanic  membrane.  To  simplify  the 
equations,  define  a  local  cartesian  coordinate  system  (x,y)  coincident  with 
the  direction  of  the  radial  and  circular  fibers.  Since  the  inner  circular 
fibers  are  approximately  orthogonal  to  the  outer  radial  fibers,  selecting  the 
local  x  axis  to  coincide  with  a  radial  fiber  will  result  in  y  running  approx¬ 
imately  along  a  circular  fiber  (Fig.  3.3). 


Figure  3.3.  Sketch  of  the  fiber  arrangement  of  the  human  tympanic  membrane. 

R:  radial  fibers;  P:  parabolic  fibers;  C:  circular  fibers  (based 
on  Shimada  and  Lim,  1971). 

As  a  first  approximation  it  is  assumed  that  the  radial  and  circular 
fibers  represent  a  majority  of  the  stiff  fibers.  The  material  is  idealized 
mechanically  as  consisting  of  a  set  of  locally  orthogonal  radial  and  circular 
fibers  imbedded  in  a  base  tissue.  Mechanical  properties  of  the  fibers  and 
base  material  are  independent  and  combine  to  give  global  stiffnesses  in  eq. 
(3.1.1).  Based  on  this  idealization,  the  bending  stiffnesses,  membrane  stiff¬ 
nesses,  and  torsional  stiffnesses  will  be  estimated  in  the  local  circumferen¬ 
tial  and  radial  directions  using  the  geometry  and  properties  of  constituent 
components. 

3.2  Membrane  Stiffness 

The  membrane  stiffness  represents  the  resistance  to  tensile  deformation 
in  the  plane  of  the  membrane.  This  stiffness  is  analogous  to  Young's  modulus 


for  isotropic  materials  multiplied  by  the  thickness.  For  the  idealized 
material.  Appendix  A  approximates  the  membrane  stiffness  in  the  direction  of 
the  i  fiber  as  (i  may  represent  only  the  radial  direction  or  the  circular 
di recti  on) 
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where  the  subscript  b  represents  the  base  material,  d  is  the  equivalent  fiber 

diameter,  h  is  the  total  thickness,  \  is  the  center  to  center  fiber  spacing, 

E  is  the  elastic  moduli  for  the  constituent,  and  v  is  Poission's  ratio  for  the 

constituent. 

Using  estimates  for  the  cat  (Funnell,  1975)  h  =  40  x  I0~4cm,  d^  = 

8  x  10-4  cm,  E(j  =  108  dyn/cm2,  E^  *  1010  dyn/cm3,  Yb  =  Yi  =  0.3,  we  find  for 
sparse  fiber  spacing  that  0-j. sparse  “  4.4  x  io4  dyn/cm.  For  tight  spacing 
with  the  fibers  touching  each  other  Qi -tight  **  6*9  *  1°6  dyn/cm.  In  the 
radial  direction,  most  of  the  membrane  is  packed  very  tightly  with  radial 
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Figure  3.4  Estimated  membrane  stiffness  vs.  radial  position. 


fibers  (Shimada  and  Lim,  1971).  In  the  circular  direction,  however,  spacing 

is  relatively  tight  near  the  annular  ring  but  becomes  very  sparse  near  the 

umbo.  Hence,  membrane  stiffness  is  highly  anisotropic  particularly  near  the 

umbo.  Figure  3.4  graphically  shows  estimates  of  the  axial  stiffnesses  in  the 

radial  and  circumferential  directions  for  various  positions  on  the  tympanic 

membrane.  If  we  were  to  calculate  the  effective  elastic  moduli  in  the  radial 

direction  near  the  umbo  based  on  the  total  thickness  the  result  would  be 

E*  =  8.6  *  109  dyn/cm2.  No  experimental  data  of  the  membrane  in  pure  axial 
m 

tension  has  been  reported,  and  hence,  we  have  no  data  to  compare  to  this 
stiffness.  The  numerical  values  used  to  calculate  the  axial  stiffness  are 
physically  consistent  with  the  types  of  materials  in  the  actual  drum.  In 
fact,  it  will  be  shown  in  the  following  section  that  the  above  values  repro¬ 
duce  elastic  moduli  derived  from  bending  stiffness  tests. 


3.3  Bending  Stiffness 

By  applying  the  Ki rchhoff  hypothesis  for  bending  of  shells,  it  is  pos¬ 
sible  to  derive  the  bending  stiffness  of  the  idealized  drum  material  in  terms 
of  its  construction  and  constituent  materials.  By  integrating  the  curvature 
induced  stress  over  the  cross  section  (see  Appendix  A)  the  bending  stiffness 
for  the  i  direction  is  estimated  to  be 
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In  addition  to  the  parameters  defined  in  (3.2.1),  a^  is  the  distance  from  the 
i  layer  to  the  neutral  surface  in  the  i  direction,  h^  is  the  thickness  above 
the  neutral  surface,  and  hj  is  the  thickness  below  the  neutral  surface. 


Fibers  of  the  eardrum  are  much  stiffer  than  the  supporting  cellular 
structure  forcing  the  neutral  surface  to  be  centered  in  the  middle  fiberous 
layer.  Hence,  for  the  eardrum  eq.  (3.3.1)  approximately  reduces  to 

0.  -  Jb  {h2  +  hi  .  £]_}  +  J1  £l_  .  (3.3.2) 
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Using  parameter  estimates  from  the  previous  section  we  find  for  a 
tightly  packed  fiber  spacing,  representative  of  the  radial  direction  near  the 
umbo,  that  the  bending  stiffness  is  approximately  1.37  dyn-cm.  This  stiffness 
remains  essentially  constant  as  the  fiber  density  is  changed  because  it 
depends  primarily  upon  the  base  material  parameters.  Insensitivity  of  the 
bending  stiffness  to  fiber  densities  explains  why  von  Bekesey  (1949),  and 
Kirikae  (1960)  experimentally  measured  a  homogeneous  bending  stiffness  value 
across  the  pars  tensa.  The  fact  that  the  bending  stiffness  is  approximately 
homogeneous  has  been  misinterpreted  by  some  researchers  as  indicating  that  the 
drum  acts  isotropically. 

Despite  the  nearly  homogeneous  bending  stiffness,  comparison  to  analy¬ 
tically  estimated  membrane  stiffnesses  indicates  that  the  pars  tensa  is  highly 
anisotropic.  To  illustrate  this,  if  the  material  is  treated  as  a  single 
phase  material  then  Young's  modulus  can  be  extracted  from  the  bending  stiff¬ 
ness  using  D.  =  E*  h2/12(l-v2).  For  Di  =  1.37  dyn-cm  and  v  =  0.3  we  find 
Young's  modulus  from  the  bending  stiffness  to  be  E^j  =  2.35  x  108  dyn/cm2.  This 
elastic  modulus  value  is  essentially  identical  to  the  moduli  experimentally 
obtained  by  several  researchers  (von  Rekesey,  1949;  Kirikae,  1970).  All  of 
the  experimental  values  were  obtained  by  measuring  the  bending  stiffness 
and  calculating  E^  from  D.  just  as  was  done  here  making  the  bending  stiffness 
estimate  of  eq.  (3.3.1)  directly  comparable  to  experimental  values.  The 


elastic  modulus  E*  obtained  from  the  membrane  stiffness,  however,  is  not 
m 

directly  comparable  to  experimentally  derived  moduli  values.  If  the  eardrum 

were  isotropic  then  Ej  would  be  identical  to  E*.  For  the  model  eardrum  ultra- 

d  m 

structure  the  membrane  based  moduli  E*  greatly  exceeds  the  bending  based 

m 

moduli  E^  (Fig.  3.5). 
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Figure  3.5.  Approximate  membrane  stiffness/bending  stiffness  vs.  total 

thickness.  Circle  indicates  TM  range,  dashed  line  represents 
isotropic  limit. 


If  the  tympanic  membrane  were  treated  as  isotropic  or  even  single  phase 
anisotropic  then  matching  the  bending  stiffness  to  experimental  values  would 
underpredict  the  membrane  stiffness.  This  fact  strongly  indicates  that  a 
single  phase  constitutive  law,  isotropic  or  orthotropic,  may  not  be  reason¬ 
able  for  the  tympanic  membrane. 


3.4  Torsional  Stiffness 

The  torsional  stiffness  may  be  estimated  using  results  of  Appendix  A 


in  a  way  similar  to  the  bending  stiffness.  From  Appendix  A 


Using  parameters  from  Section  3.2  the  torsional  stiffness  is  approximately 
S12  ~  0.19  dyn  -  cm.  If  the  material  were  single  phase  then  Young's  moduli 

determined  by  E*  h3/12(l-v2)  ~  (D^v^  +  S12^  wou1d  be  Es  ~  1,3  x  108  dyn/cm2* 
This  modulus  differs  from  the  bending  derived  modulus  by  a  factor  of  2,  and 

differs  from  the  membrane  derived  modulus  by  an  order  of  magnitude.  This, 
again,  discredits  use  of  the  single  phase  constitutive  laws  for  the  eardrum. 

3.5  Conclusions 

In  addition  to  membrane,  bending,  and  torsional  terms  some  coupling 
between  these  terms  may  exist  in  the  tympanic  membrane.  Selecting  a  locally 
orthogonal  coordinate  system  approximately  coincident  with  the  fiber  direc¬ 
tions  eliminates  most  of  the  coupling.  Since  the  fibers  are  not  exactly 
orthogonal  and  since  the  radial  and  circular  fibers  are  on  different  levels, 
some  coupling  will  exist.  The  ultrastructure,  however,  suggests  that  the 
coupling  is  small  and  that  it  may  not  be  necessary  to  include  such  terms  to 
model  primary  action  of  the  drum.  For  this  reason,  the  coupling  moduli  will 
not  be  estimated  here. 

Example  stiffness  calculations  are  done  in  Sections  3.2,  3.3  and  3.4 
for  the  pars  tensa  of  a  cat.  By  using  different  parameters,  the  equations  may 
also  be  applied  to  the  pars  flaccida.  The  stiffnesses  are  local  and  change 
depending  on  the  local  ultrastructure  of  the  tympanic  membrane.  By  changing 
parameters  the  model  ultrastructure  is  changed  geometrically  and  materially. 
Parameter  selection  allows  the  stiffness  formul i  to  estimate  mechanical  pro¬ 
perties  for  other  mammals  including  guinea  pigs,  sheep,  and  humans. 


In  conclusion  it  is  stressed  that  the  fiberous  tympanic  membrane  ultra¬ 
structure  has  anisotropic  mechanical  properties.  Evidence  of  isotropic 
bending  stiffness  should  not  be  misconstrued  as  evidence  of  isotropy.  In 
fact,  estimation  of  membrane  and  torsional  stiffnesses  indicates  very  aniso¬ 
tropic  properties.  The  anisotropy  is  primarily  seen  in  membrane  stiffness 
terms.  Modeling  the  tympanic  membrane  as  isotropic,  or  single  phase  anisotro¬ 
pic  may  induce  order  of  magnitude  errors  in  the  membrane  terms  even  if  the 
bending  stiffness  is  matched  exactly.  The  stiffnesses  derived  herein,  may 
aid  in  development  of  a  consistent  tympanic  membrane  model  by  allowing  direct 
nondimensionalization  of  anisotropic  curvilinear  shell  equations  accompanied 
with  direct  deduction  of  primary  terms. 
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3.6  Appendix  A:  Approximate  Stiffness  Decomposition 

Bending,  membrane,  and  torsional  stiffnesses  of  the  tympanic  membrane 
can  be  estimated  based  on  mechanical  properties  of  the  lamina  propri,  mucous, 
and  epidermal  layers.  To  derive  approximate  expressions  for  stiffnesses  in 
terms  of  constituent  components,  consider  an  idealized  ultrastructure  consist¬ 
ing  of  locally  orthogonal  fibers  imbedded  in  a  base  material  as  shown  in 
Figure  3.6.  The  composite  shell  is  constructed  of  a  set  of  fibers  running 
parallel  to  the  x  axis  having  diameters  d^  and  spacing  Xj,  and  a  set  of  fibers 
running  parallel  to  the  y  axis  having  diameters  d£  and  spacing  X2.  The 
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Figure  3.6  Local  orthogonal  fiber  composite  model  of  the  tympanic  membrane 
material  (Fiber  spacing  and  size  not  to  scale.) 
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relatively  stiff  orthogonal  set  of  fibers  is  imbedded  in  a  base  material 
having  thickness  h.  The  force  and  moment  equilibrium  equations  must  hold 
regardless  of  the  presence  of  the  fiber  reinforcement  in  the  drum.  Hence,  in 
order  to  formulate  continuous  shell  equations  that  approximate  behavior  of  the 
discrete  composite  structure,  we  need  to  determine  material  constants  or 
stiffnesses  that  will  represent  the  discrete  system.  One  way  to  estimate  the 
stiffnesses  would  be  to  determine  the  behavior  of  one  unit  element  of  the  grid 
and  substitute  the  resulting  parameters  into  equilibrium  equations.  (Fig. 
3.7).  Renton  (1970),  Lightfoot  (1964)  and  Yettram  (1964)  have  applied  this 


Figure  3.7  Equilibrium  of  a  grid  element.  Njj  denotes  forces  and  M-jj 
denotes  moments  on  each  cross  sectional  face. 


approach  to  estimate  parameters  of  plates  for  finite  difference  gridworks. 
When  discretizing  a  continuous  plate  into  a  finite  gridwork  particular  atten¬ 
tion  should  be  paid  to  the  effect  of  Poisson's  ratio.  Obviously,  if  material 
coupling  between  the  x  and  y  directions  is  removed  then  increased  stiffness 
due  to  Poisson's  ratio  will  be  lost.  This  fact  can  bseen  by  inspecting 
cross  product  terms  in  the  orthotropic  plate  equations.  For  a  homogeneous 
orthotropic  plate,  the  coefficient  of  the  product  moment  term  (representing 
torsional  resistence)  is 

Dx  vxy  +  Dy  Vyx  +  Sxy  +  SyX  , 

where  Dx  is  the  bending  stiffness  in  the  x  direction,  Dy  is  the  bending  stiff¬ 
ness  in  the  y  direction,  vxy  is  the  first  Poisson  ratio,  vyx  is  the  second 
Poisson  ratio,  and  SXy  are  torsional  stiffnesses. 

For  a  discrete  gridwork  the  coefficient  multiplied  by  Poisson's  ratio 
will  approach  zero  as  dn/xn  approaches  zero.  Hence,  if  we  have  a  plate  with 
uniform  bending  stiffness  D,  and  we  have  a  gridwork  with  the  same  bending 
stiffness  D  then  the  overall  torsional  stiffness  of  the  plate  will  exceed  that 
of  the  gridwork.  The  same  comments  apply  to  the  effect  of  Poisson's  ratio  on 
axial  stiffnesses,  such  that  the  deflection  of  a  gridwork  having  axial  stiff¬ 
ness  and  bending  stiffness  identical  to  that  of  the  plate,  may  deflect  quite 
differently  from  the  plate  under  identical  loads.  With  these  comments  in  mind 
we  would  like  to  find  the  appropriate  stiffnesses  for  use  in  an  equivalent 
continua  representation  of  the  discrete  composite  plate. 

In  order  to  estimate  the  bending  stiffness  in  the  x  direction,  suppose 
we  assume  that  the  Kirchhoff  hypothesis  is  reasonable,  such  that  the  bending 


moment  may  be  obtained  as  an  integral  of  the  stress  over  the  section. 


mX2  h2 

M  =  Lim  __L  /  f  a  zdz  dy 
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Assuming  linear  Kirchhoff  type  strain,  the  integral  may  be  estimated  in  the 


f ol  1  owi  ng  way 
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where  w  is  transverse  deformation  and  x,  y,  z  subscripts  on  w  denote  partial 
differentiation.  It  has  been  assumed  that  each  of  the  constituent  components 
behave  in  a  linear  elastic  manner.  The  subscripts  appearing  on  the  material 
constants  indicate:  1-fibers  in  the  first  direction,  2-fibers  in  the  second 
direction,  b-base  material.  Figure  3.B  indicates  definitions  of  thickness  and 
spacing  parameters. 

In  addition  to  the  above  form  the  bending  moment  can  be  written  as 
M  =  °2^wxx  +  v2i  wyy^  de^in’n9  the  bending  stiffness  02  and  Poisson's  ratio  as 
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Performing  the  integration  for  Kirchhoff  type  strain  gives 


(3.6.6) 
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Once  again  an  interchange  of  subscripts  gives  S12. 

The  axial  membrane  stiffness  may  be  derived  by  integrating  the 
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membrane 


stress  over  the  thickness  of  the  material.  In  the  x  direction  the  in-plane 


membrane  force  is 
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where  u  is  deformation  in  the  x  direction,  v  is  deformation  in  the  y  direction 
and  is  the  membrane  stiffness.  Upon  integration 
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Finally  the  in-plane  shear  is  estimated  by  integrating  the  shear  stress 


over  the  thickness 


(3.6.10) 
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Performing  the  integration  we  find 
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Interchanging  subscripts  1  and  2  will  provide  Q12* 

With  these  stiffnesses  the  equivalent  conti nua  problem  has  been  formu¬ 
lated.  All  of  the  expressions  are  approximate  and  are  valid  only  for  a  limited 
range  of  parameters  and  conditions.  When  applied  to  dynamic  problems,  the 
expressions  are  reasonable  only  for  low  frequency  vibrations  where  material 
wave  propagation  is  negligible.  An  alternative  to  the  "strength  of  materials" 
approach  presented  herein  is  presented  by  Tiersten  and  Jahanmir  (1977)  modeling 
composites  as  interpenetrating  conti nua. 


CHAPTER  4 


SIMPLE  EVIDENCE  OF  MEMBRANE  STRESS  AND  DAMPING  IN  THE  TYMPANIC  MEMBRANE 
4.1  Introduction 

In  order  to  formulate  a  simple  yet  physically  realistic  model  of  the 
tympanic  membrane  we  must  study  its  structure  and  behavior  in  order  to  identify 
the  important  features.  The  problem  is  to  properly  identify  all  of  the 
dominant  mechanisms.  There  are  two  basic  schemes  of  searching  for  important  or 
primary  features  when  formulating  a  model.  The  first  method  is  to  write  a 
complete,  and  relatively  large,  set  of  equations  describing  behavior  of  the 
eardrum  and  subsequently  use  known  behavior  to  nondimensionalize  and  simplify 
the  model.  For  the  eardrum,  this  would  involve  writing  von  Karman  type  shell 
equations  for  the  fiber  composite  structure  in  general  curvilinear  coordinates 
coupled  with  acoustic  models  of  connecting  fluids  and  ossicular  chain  constraints. 
The  curvilinear  geometry,  fiber  composite  structure,  and  internal  damping  of  the 
tympanic  membrane  would  make  even  a  numerical  solution  impractical  at  this  time. 
Although  it  is  possible  to  formulate  a  relatively  complete  model,  many  features 
may  be  of  secondary  importance.  If  such  terms  could  be  located  apriori,  then 
the  modeling  task  could  be  greatly  simplified.  While  the  first  method  of 
searching  for  primary  terms  requires  a  relatively  complete  formulation,  the 
second  method  attempts  to  exclude  selected  physical  features  by  estimating 
their  contribution  prior  to  formulation  of  the  model.  In  the  interest  of  sim¬ 
plicity,  importance  of  several  physical  features  of  the  tympanic  membrane  will 
be  estimated  for  consideration  in  an  eardrum  model.  Actual  experimental  results 
will  be  used  to  make  the  analytical  estimates. 

Since  the  tympanic  membrane  is  a  thin  fiber  composite  shell  it  is  natural 
to  think  of  it  in  terms  of  dynamical  shell  theory.  When  formulating  a  shell 
model,  the  physics  may  be  divided  into  three  catagories.  1)  Inertia:  transverse. 


rotary,  shear,  axial  ...  2)  Restoring  forces:  shear,  bending,  axial,  external 

...  3)  Damping:  transverse,  bending,  shear  ...  .  This  chapter  addresses 

existence  of  dominant  mechanisms  in  each  of  these  catagories  beginning  with  a 
discussion  of  inertia. 

Using  holographic  interferometry  the  vibrational  shape  of  the  tympanic 
membrane  has  been  recorded  under  various  experimental  conditions  (Khanna  and 
Tornndorf,  1972;  Tonndorf  and  Khanna,  1972;  Dancer,  1975;  Ogura,  1974).  The 
displaced  shapes  can  be  used  to  estimate  magnitudes  of  shear,  transverse,  and 
rotary  inerta  mechanisms.  Most  of  the  membrane  motion  is  transverse  to  its 
resting  position  indicating  that  transverse  inertia  dominates  other  inertia 
types  (for  frequencies  in  the  audible  range).  The  only  possible  exception  is 
in  the  area  of  taught  radial  fibers  separating  the  anterior  from  posterior 
regions  where  connection  of  the  malleus  induces  larger  than  average  tangential 
accelerations.  However,  since  the  audio  excitation  frequencies  represent 
speeds  much  slower  than  the  sound  speeds  of  shear  and  dilatational  waves, 
local  material  accelerations  represent  much  less  kinetic  energy  change  than 
global  accelerations  (Graff,  1970).  Therefore,  rectilinear  momentum,  primar¬ 
ily  tranverse,  represents  the  primary  inertia  in  the  vibrating  tympanic  mem¬ 
brane. 

Determination  of  the  dominant  restoring  mechanisms  is  not  as  simple  as 
the  inertia  mechanism.  Based  on  some  kinematic  considerations,  however,  it  is 
possible  to  estimate  the  strain  energy  density  stored  by  bending,  the  potential 
energy  density  stored  by  membrane  stress,  and  the  potential  energy  stored  by 
shear  mechanisms.  In  order  to  estimate  the  energy  density  stored  in  membrane 
type  action  it  is  necessary  to  first  determine  the  magnitude  of  an  in-plane 
stress  existing  in  the  eardrum.  There  are  several  ways  to  infer  the  magnitude 
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of  the  membrane  stress  based  on  experimental  results  -  three  methods  will  be 

presented  here  and  compared. 

4.2  Approximate  Global  Static  Equilibrium 

Consider  an  idealized  boundary  that  supports  the  tympanic  membrane 

through  membrane  type  forces  as  indicated  in  Fig.  4.1.  In  reality,  there  will 

exist  some  transverse  shear  and  bending  moment  at  the  boundary  accompaning  the 

membrane  stress,  however  the  intent  here  is  to  estimate  the  membrane  stress 

if  it  dominates  over  the  bending  stresses.  To  obtain  an  order  of  magnitude 

approximation,  consider  a  cone  of  diameter  2a,  and  apex  angle  2a  forced  by  a 

pressure  P,  tympani  force  FT,  ossicular  moment  M  and  membrane  stress  a*. 

'i  m  m 

Summing  moments  about  the  manubrium  instant  center  of  velocity  we  find 

q*  ~  Ftb  +  P*a3  -  Mm  ,  (4.2.1) 

m  0.17  a*h  sin  a 


(a,0) 


Figure  4.1  Approximate  global  equilibrium  of  the  tympanic  membrane. 
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where  the  moment  Mm  may  include  reverse  inertial  terms  representing  an  inte¬ 
grated  mass  of  the  hearing  system  as  well  as  damping  and  compliance  terms. 
Also,  writing  the  moment  as  an  operator  on  malleus  displacement  then 
Mm  =  <(> K  [<j> , <(> where  $  is  the  instantaneous  rotation  of  the  malleus. 

The  tensor  tympani  force  and  perhaps  some  of  the  moment  Mm  can  be 
viewed  as  inducing  a  prestress  ap  such  that  a*  is  the  sum  of  a  prestress  and 
a  deformation  induced  stress.  To  estimate  the  magnitude  of  the  membrane 
stress  a*,  dimensions  representative  of  a  cat  are  combined  with  observed  de¬ 


formation  to  obtain  the  following 


~  on  +  1.23  x  in5  dyn  . 

r  cm* 


(4.2.2) 


The  estimate  in  (4.2.2)  was  done  at  100  db  SPL  using  deformation  rotation 

<t>  ~  1.5  x  io"5  rad.  measured  by  Khanna  and  Tonndorf  (1972)  and  h  =  40  x  io-4  cm, 

a  =  0.45  cm,  b  =  0.25  cm,  a  =  0.64  rad,  and  from  Funnel  1  (1977)  K  -  28  x  103 

m 

dyn  -  cm.  The  membrane  stress  estimate  of  (4.2.2)  is  insensitive  to  any  errors 
in  the  malleus  spring  constant.  Also,  the  prestress  ap  induced  by  the  tensor 
tympani  or  tendons  has  been  left  arbitrary  for  now. 

4.3  Radial  Fiber  Strain 

To  verify  the  estimate  of  radial  stress  based  on  global  static  equili¬ 
brium  (4.2.2),  strain  in  the  fibers  can  be  estimated  from  experimental  results 
and  combined  with  the  membrane  elastic  moduli  to  obtain  a  comparative  value. 

The  strain  of  a  radial  fiber  for  small  deformations  can  be  written  as  a  sum  of 
strain  due  to  endpoint  motion  and  strain  induced  by  transverse  deformation 
er  “  where  the  superscript  e  represents  endpoint  motion  and  the 

superscript  q  represents  quasilinear  stress  induced  by  the  transverse  deforma¬ 
tion. 


To  estimate  strain  due  to  endpoint  motion  consider  boundaries  defined 


by  a  cone  shape  such  that  any  planar  cut  through  the  umbo  perpendicular  to  the 
plane  of  the  annular  ring  will  yield  a  triangular  cross  section.  The  cross 
section  is  allowed  to  change  with  e  so  the  cone  may  be  distorted.  This  shape 
is  used  to  estimate  endpoint  motion  only,  so  curvature  of  the  drum  does  not 
come  into  the  calculation.  Fig.  4.2  illustrates  the  approximate  boundaries 
under  consideration. 


Figure  4.2  Membrane  boundaries  used  to  estimate  radial  fiber  strain. 


Strain  on  a  radial  fiber  due  to  endpoint  motion  is  given  by 


2  2 
e(e)  _  1  Cs  _  C-j 

r  " 
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(4.3.1) 


where  C..  and  C$  are  initial  and  final  fiber  lengths>  respectively.  Noting 
that  b  and  l  are  constant  during  the  deformation,  and  $  is  small  gives 


where  $  and  aQ  are  projections  onto  the  plane  of  the  triangle  containing  the 
radial  fiber  of  interest,  while  l,  aQand  b  are  shown  in  Fig.  4.1.  Also 
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We  must  add  to  (4.3.2)  the  strain  due  to  transverse  deformation.  Since  the 
drum  is  initially  convex,  strain  caused  by  positive  endpoint  motion  will  be 
positive,  and  strain  due  to  inward  transverse  deformation  will  be  negative. 
If  we  denote  u(r)  as  the  transverse  deformation  then 
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dr 


(4.3.5) 


Approximating  u(r)  for  low  frequency  vibration  (less  than  1000  Hz)  as 
u  ~  %  sin  0(r  -  si )/(s2  “  si))  9ives  ej.q)  -  6.9[u0/(S2  -  Sj]2,  where  uq  is 
the  maximum  transverse  deflection  of  the  drum  relative  to  a  cone  deformed  to 
meet  the  endpoint  deformation. 

In  addition  to  the  values  representing  the  cat  geometry  given  in  Sec¬ 
tion  4.2,  at  100  dB  SPL  the  maximum  transverse  displacement  after  correcting 
for  endpoint  motion  is  approximately  uQ  ~  0.21  *  10-4  cm.  Using  typical 
values  we  find  that  the  strain  associated  with  endpoint  motion  is  approximat¬ 
ely  -  1.17  x  io-5  and  the  strain  induced  by  transverse  deformation  is 
approximately  ~  -2.3  x  1Q-8.  Hence,  the  compressive  strain  induced  by 


transverse  deformation  and  change  in  curvature  is  two  orders  of  magnitude  less 
than  the  tensile  strain  induced  by  motion  of  the  endpoints.  This  observation 
discredits  simple  Helmholtz  curved  lever  theories  of  membrane  action.  The 
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curved  lever  theories  require  existence  of  a  compressive  stress  in  the  radial 

fibers  under  the  conditions  estimated  here.  The  fact  that  a  tensile  strain 

exists,  not  a  compressive  strain,  suggests  the  curved  lever  mechanism  is  a 

secondary  effect  and  does  not  dictate  drum  function. 

To  verify  the  fiber  strain  estimate,  the  membrane  stiffness  estimated 

by  equation  (3.2.1)  can  be  used  to  estimate  the  stress  a**.  Doing  so,  we  find 

a**  ~  1.1  x  105  dyn/cm2  +  cr_  .  This  result  is  essentially  identical  to  that 
m  P 

obtained  by  global  equilibrium  (4.2.2)  and  hence  seems  very  reasonable. 

Since  the  membrane  stress  appears  as  an  integral  of  the  deformation, 
plus  a  prestress,  it  will  be  a  quasi  linear  term  in  the  tympanic  membrane  oper¬ 
ator.  If  the  prestress  op  greatly  exceeds  the  quasi  linear  term,  then  the 
system  equations  could  be  linearized  (in  terms  of  the  membrane  type  stress) 
around  the  prestress.  As  discussed  later,  some  dynamical  and  static  experi¬ 
mental  results  indicate  that  the  eardrum  behaves  as  if  a  prestress  dominates. 

4.4  Restoring  Potential  Energy  Densities 

Knowing  an  approximate  value  of  the  membrane  stress  in  typical  radial 
fibers  allows  us  to  estimate  the  potential  energy  density  stored  in  membrane 
type  deformation.  Again,  using  the  cat  as  an  example,  the  radial  fiber  along 
which  the  maximum  deflection  occurs  at  low  vibration  frequencies  has  a  shape 
similar  to  a  sine  bay.  As  a  rough  estimate,  transverse  deformation  from  the 
resting  position  for  a  radial  fiber  is  written  u(r,t)  ~  uo(t)sin(ir(r-S1)/(S2-S1)). 
Using  this,  the  average  potential  energy  per  unit  surface  area  (Um)  stored  along 

the  radial  fibers  due  to  membrane  type  forces  is  approximately 
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At  100  dB  SPt  using  estimates  for  the  cat  and  stress  a*  of  Section  4.2  we  find 
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U  -  3.02  x  10-6  dyn-cm  +  n 
m  cmz  ">P 


(4.4.2) 


where  U  _  is  due  to  prestress  in  the  tympanic  membrane  and  U  -  U  n  is  induced 
mP  m  mP 

by  deformation  alone. 

The  strain  energy  stored  in  bending  may  be  calculated  for  the  same  vibra¬ 
tional  shape  and  amplitude  as 
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Thus,  using  identical  values  for  a  cat  at  100  dB  SPL 


UD  ~  4.51  x  io-7  dyn-cm  . 

B  CTTV2 


(4.4.4) 


Based  on  these  estimates  at  100  dB  SPL,  the  membrane  type  potential  energy  is  at 
least  an  order  of  magnitude  larger  than  the  bending  strain  energy.  If  a  substan¬ 
tial  prestress  exists  in  the  eardrum,  or  is  imposed  externally,  then  the  membrane 
type  term  may  be  several  orders  of  magnitude  larger  than  the  bending  term.  The 
ratio  of  the  two  terms  provides  a  measure  of  their  relative  importance  and  is 
defined  by  2 
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(4.4.5) 


The  above  comparison  was  applied  to  bending  in  the  radial  direction  only; 
the  same  type  of  analysis  may  be  applied  in  the  circumferential  direction.  In 
order  to  estimate  potential  energy  density  due  to  membrane  action  in  the 
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circumferential  direction,  it  is  necessary  to  estimate  stress  in  the  circumferen¬ 
tial  direction.  This  can  be  done  by  considering  local  dynamic  force  balance  of 
the  membrane.  Local  principle  radii  of  curvature  plays  a  primary  role  in  relating 
radial  stress  to  circumferential  stress.  Upon  estimating  bending  and  membrane 
energy  densities  stored  due  to  deformation  in  the  circumferential  directions,  it 
is  found  that  both  are  small  in  comparison  to  membrane  energy  in  the  radial 
di rection. 

These  simple  energy  density  calculations  are  based  on  holographic 

interferometry  measurements  made  by  Khanna  and  Tonndorf  (1972),  and  hence  are 

based  on  accepted  behavior  of  the  drum.  The  only  analytical  approximation  that 

may  contain  any  appreciable  error  is  the  estimated  quasi  linear  radial  stress  o  . 

m 

The  radial  (or  meridian)  stress  used  was  estimated  using  two  distinct  methods  re¬ 
sulting  in  the  same  value  of  stress.  Both  methods  determined  the  component  of 
meridian  stress  superimposed  on  any  resting  tension  or  prestress.  A  numerical 
value  of  prestress  could  not  be  obtained  by  either  method.  If  a  prestress  exists 
then  the  radial  membrane  type  restoring  mechanism  may  dominate  by  several  orders 
of  magnitude  rather  than  just  one  order  of  magnitude. 

4.5  Simple  String  Analogy 

The  question  of  prestress  can  be  partially  resolved  by  study  of  the  vibra¬ 
tional  shape  of  the  membrane  along  with  motion  of  the  malleus.  To  address  this 
question,  suppose  we  look  at  behavior  of  a  very  simple  string  subject  to  a  radial 

stress  o  ,  supported  by  a  spring  on  one  end,  and  fixed  at  the  other  end.  This 
m 

string,  in  a  very  crude  way,  represents  a  bundle  of  radial  fibers  connected  to  a 
moving  manubrium  and  fixed  at  the  annular  ring.  (See  Fig.  4.9,  Appendix  A).  The 
string  analogy  is  by  no  means  complete  but  will  suffice  to  illustrate  several 


The  problem  of  interest  is  to  allow  the  forcing  fequency  to  vary  and 


observe  response  of  the  string  (analogous  to  TM)  and  associated  motion  of  the 

spring  (analogous  to  malleus  motion).  Solution  to  the  linear  string  problem  with 

a  spring  constraint  can  be  solved  exactly.  For  the  case  of  no  damping,  response 

of  the  string  is  shown  in  Figs.  4.3  and  4.4.  In  terms  of  infering  tympanic 

membrane  radial  prestress,  peaks  on  the  MAX  TM  AMPLITUDE/MALLEUS  AMPLITUDE  VS. 

FREQUENCY  graph  correspond  to  even  numbered  natural  frequencies  of  vibration. 

Valleys  on  the  graph  correspond  to  odd  natural  modes  of  vibration.  Since  the 

natural  frequencies  are  directly  related  to  the  stress  cr  ,  measurement  of  the 

m 

peaks  allows  calculation  of  the  radial  tensile  stress  a  . 

m 


Figure  4.3  Undamped  vibrations  of  a  string  with  one  end  supported 
by  a  transverse  spring  for  the  case  of  no  damping. 


Figure  4.4  String  analogy.  Endpoint  amplitude  /  maximum  amplitude 
for  no  damping  vs.  frequency. 

Depending  on  the  frequency,  the  point  of  maximum  displacement  on  the 
tympanic  membrane  moves.  To  illustrate  this,  at  1000  Hz  the  string  has  a  maximum 
displacement  approximately  half  way  between  the  fixed  end  (annular  ring)  and 
moving  spring  (malleus).  At  4000  Hz  the  maximum  amplitude  position  moves  toward 
the  spring.  Very  similar  behavior  is  seen  when  measuring  vibratory  behavior  of 
real  tympanic  membranes.  Manley  and  Johnstone  (1974)  have  recorded  displacement 
of  the  tympanic  membrane  of  guinea  pigs  for  several  fixed  points  over  a  wide 
range  of  frequencies  illustrating  this  behavior  (see  Fig.  4.5). 
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Figure  4.5  Displacement  of  sites  E,  F,  and  G  (see  inset  diagram  for  eardrum 
measurement  locations)  relative  to  site  D  (in  decibels)  as  a 
function  of  frequency.  The  zero  line  represents  the  displacement 
at  D  at  all  frequencies.  Relative  displacement  of  E  shown  as 
continuous  line,  filled  circles;  of  F  as  dotted  line,  open 
triangles;  and  of  G  as  dashed  line,  filled  squares.  Data  for  E 
from  one  individual,  for  F  is  a  mean  of  three  individuals,  and 
for  G  is  a  mean  of  two  individuals.  The  trends  in  all  cases  were 
the  same.  From  Manley  and  Johnstone  (1974). 


By  careful  inspection  of  Manley  and  Johnstone's  data,  keeping  in  mind 

f 

that  the  position  of  the  maximum  displacement  moves  with  frequency,  it  is 
possible  to  locate  the  first  four  peaks  on  the  TM/Malleus  amplitude  graph. 
They  appear  approximately  at  a>  ~  1.86  KHZ,  u  -  3.9  KHZ,  a>  -  5.1  KHZ,  and 


a)  -  7.4  KHZ.  These  peaks  correspond  roughly  to  even  modes  of  vibration  of  the 
8 

radial  fiber  on  which  data  was  taken.  If  we  match  the  value  of  u>2  measured  for 

s 

the  guinea  pig  to  the  result  for  the  simple  string  problem  we  find  u>2  ~  1.86  KHZ 
s  s  s 

<D  ~  3.7  KHZ,  u»  -  5.6  KHZ,  u>  ~  7.4  KHZ.  The  superscript  s  indicates  the  string 
*+6  8 

analogy.  This  agreement  is  outstanding.  Any  differences  are  smaller  than 
experimental  accuracy  of  the  measured  frequencies.  Similar  results  can  be 
obtained  for  the  cat  eardrum.  A  graph  for  the  cat  by  Khanna  and  Tonndorf  (1972) 
is  shown  in  Fig.  4.6. 


Figure  4.6  Maximum  TM  amplitude  vs.  malleus  amplitude. 
(From  Khanna  and  Tonndorf,  1972) 


By  matching  peaks  on  the  TM/Malleus  amplitude  experimental  graphs  to  the 
string  problem  it  is  found  that  the  radial  stress  must  be  of  order  106  dyn/cm2 


\  vV.’ 


for  the  cat.  As  an  indirect  confirmation  of  this  stress.  Funnel  1  (1975)  used  a 


prestress  of  4  x  106  dyn/cm2  in  a  single  phase  plane  membrane  model  to  produce 
complex  vibratory  shapes  observed  in  cats.  This  value  of  stress  is  an  order  of 
magnitude  larger  than  the  quasi  linear  stress  previously  obtained  in  (4.2.2) 


indicating  that  there  exists  a  prestress  in  the  cat  eardrum  of  order  106  dyn/cm2. 
Since  the  tensor  tympani  plays  a  role  in  determining  the  prestress,  the  value  may 
be  sound  pressure  level  dependent.  The  106  dyn/cm2  is  based  on  results  at  100  dB 


SPL. 

In  addition  to  the  location  of  the  peaks,  a  second  item  of  importance  is 
spacing  of  the  natural  frequencies.  By  experimental  observation  it  is  seen  that 
peaks  on  the  TM/MALLE1JS  amplitude  graph  have  spacing  of  "2n"  where  n  =  1,2,3,..., 
corresponding  to  a  natural  frequency  have  spacing  of  "n" .  This  spacing  can  be 
reproduced  by  a  mathematical  model  if  the  model  contains  a  dominant  membrane  type 
restoring  mechanism.  If  the  mathematical  model  contains  only  bending  restoring 
mechanisms,  then  the  spacing  of  the  peaks  would  be  approximately  (2n)2.  Spacing 
based  on  n2  is  not  seen  experimentally  discrediting  plate  or  shell  'bending 
models'  of  the  tympanic  membrane. 


4.6  Perforation  Tests 

As  additional  evidence  of  existence  of  a  tension  field  in  the  tympanic 
membrane,  perforation  tests  are  often  mentioned  (Kirikae,  1960;  Funnel  1  ,  1975). 

A  perforation  test  is  done  by  cutting  the  membrane  and  measuring  post  cut  defor¬ 
mation.  Due  to  tissue  damage  caused  by  the  cut  along  with  post  cut  tissue 
distortion,  it  is  questionable  to  attempt  to  derive  a  numerical  value  of  stress 
from  such  observations.  To  address  this  question  suppose  we  assume  for  the 
moment  that  no  tissue  distortion  exists  and  apply  a  "no  distortion"  model  to  the 
slit  test  problem. 

Modeling  the  slit  as  a  Mode  1  crack,  deformation  (v)  perpendicular  to  the 
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where  o0  is  the  prestress,  <j>  is  the  angular  coordinate;  r  is  the  radial  coor¬ 
dinate,  G  is  the  modulus  of  rigidity,  E  is  Young's  modulus,  v  is  Poisson's 
ratio,  and  a  is  the  crack  length.  Using  material  parameters  derived  in 
Chapter  3  for  a  fiber  composite  cat  eardrum  drum,  and  slit  test  results  by 
Kirikae  (1960)  we  have  v  -  0.05  mm,  a  -  0.25  mm,  E  ~  8.6  x  109  dyn/cm2.  These 
values  yield  a  prestress  of  o0  -  4.9  x  107  dyn/cm2.  This  value  is  much  larger 
than  the  value  inferred  from  vibrational  behavior  of  the  drum  and  is  believed 
to  be  in  error  due  to  tissue  distortion.  Kirikae  observed  that  the  slits 
opened  slowly  over  a  time  of  approximately  ten  seconds.  The  initial  deforma¬ 
tion  may  have  indeed  been  elastic  and  due  to  membrane  type  stress.  The 
subsequent  slow  deformation  due  to  tissue  degredation  causes  calculations, 
such  as  the  one  above,  to  be  in  error.  More  detailed  real  time  records  of 
post  cut  tissue  deformation  may  allow  future  researchers  to  obtain  stress 
values  from  slit  test  data.  The  data  presently  available,  however,  can  only 
be  used  to  make  qualitative  comparisons  and  indicate  higher  stress  in  the 
radial  fibers  than  in  the  circular  fibers  (as  expected  from  calculations  in 
Sections  4.2  and  4.3) . 

4.7  Dissipation  and  Damping 

Evidence  of  damping  within  the  structure  of  the  tympanic  membrane  is 
indicated  by  several  experimental  results.  Moller  (1965)  recorded  a  phase 
shift  in  average  displacement  of  the  eardrum  relative  to  the  stimuli  for 
anesthetized  cats  having  immobilized  ossicular  chains.  If  the  ossicular  chain 
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is  fixed  then  the  phase  shift  seen  in  acoustic  impedance  measurements  is  inde¬ 
pendent  of  the  cochlea  and  ossicular  chain.  Hence,  damping  indicated  by  the 
phase  shift  must  be  in  the  tympanic  membrane,  fluid  of  the  tympanic  cavity, 
or  air  in  the  outer  ear.  The  absence  of  appreciable  flow  within  the  outer 
ear  and  tympanic  cavity  suggests  that  structural  damping  within  the  tympanic 
membrane  plays  a  primary  role. 

In  addition  to  the  phase  shift  data,  the  vibrational  amplitude  of  the 
tympanic  membrane  can  be  used  to  infer  internal  damping  within  the  eardrum. 
Accompanying  the  phase  shift  Moller  (1965)  observed  a  rapid  decrease  in 
acoustic  impedance  for  cats  with  immobilized  ossicular  chains.  This  impe¬ 
dance  decrease  infers  an  average  amplitude  decrease  that  would  not  be  as  pro¬ 
nounced  without  damping. 

Supporting  the  global  impedance  implications,  Manley  and  Johnston  (1974) 
directly  measured  displacements  of  various  points  on  the  tympanic  membrane  of 
guinea  pigs  using  the  Mossbauer  technique.  Their  results  show  a  very  clear 
decrease  in  malleus  amplitude  over  the  entire  audible  frequency  range.  For 
frequencies  up  to  20  KHZ  they  also  show  that  the  ratio  of  peak  tympanic  mem¬ 
brane  displacement  to  malleus  displacement  does  not  display  the  same  decrease 
and  in  fact  shows  a  discrete  number  of  increasing  peaks.  The  decrease  in 

I 

malleus  displacement  with  increasing  frequency  that  is  accompanied  by  a  sub¬ 
stantial  decrease  in  tympanic  membrane  vibration  amplitude  directly  indicates 

existence  of  internal  structural  damping  within  the  eardrum. 

I 

To  include  damping  and  dissipation  in  the  model  of  the  tympanic  membrane 
it  is  necessary  to  determine  the  type  of  damping.  Most  continuum  plate  or 
|  shell  models  that  include  damping  use  a  "transverse  type"  mechanism  having  a 

dissipative  force  proportional  to  the  transverse  velocity.  This  type  of 
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damping  however,  has  no  direct  physical  connection  to  the  structure  of  the 
eardrum.  It  may  be  possible  to  crudely  equate  the  energy  dissipation  caused 
by  the  transverse  velocity  term  to  that  observed  in  the  eardrum,  but  such  an 
endeavor  violates  the  physics  of  the  problem  by  introducing  nonphysical  terms 
i n  the  model . 

Rather  than  blind  use  of  a  transverse  velocity  type  dissipation,  the 
structure  of  the  eardrum  suggests  a  "bending  type"  dissipation.  To  show  how 
the  tympanic  membrane  structure  gives  rise  to  a  bending  type  damping,  suppose 
we  look  at  behavior  of  a  three  layer  lamina  representing  the  membrane  (Lim, 
1968a,  1968b,  1970). 


Undsfonad  state 


Dcforaad  state 


Neutral  surface 


Epidermal  layer 


Lamina  propria 


Mucous  layer 


Figure  4.7  Bending  dissipation  in  the  epidermal  and  mucous  layers. 

Bending  of  the  tympanic  membrane  causes  the  mucous  and  epidermal  layers 
to  change  shape  (maintaining  similar  volume  as  indicated  in  Fig.  4.7).  Since 
the  mucous  and  epidermal  layers  are  water  intensive  cellular  structures,  the 
force  required  to  change  their  shape  will  increase  with  the  speed  of  displace¬ 
ment.  By  looking  at  changes  in  curvature  induced  by  deformation  it  is  found 


that  linear  bending  type  dissipation  for  small  deformations  is  proportional  to 
two  spatial  and  one  time  derivative  of  the  deformation  (denoted  uxxt).  The 
"bending  type"  damping  (uxxt)  differs  from  transverse  type"  (ut)  damping 
essentially  by  changing  the  slope  on  the  amplitude  vs.  frequency  graph. 

To  illustrate  the  difference  consider  two  strings,  one  with  transverse 
damping  and  one  with  bending  type  damping.  Asymptotic  behavior  of  the  strings 
can  be  found  using  a  multi  scale  method.  In  order  to  compare  the  results  note 
that  Fig.  4.8  has  a  log  frequency  scale,  so  equally  spaced  natural  frequencies 
will  become  grouped  together  in  the  high  frequency  range  indicating  that 
forcing  near  a  natural  frequency  is  appropriate  for  comparison  to  the  graph  at 
high  frequencies.  At  low  frequencies  the  eigenvalues  are  separated  by  sub¬ 
stantial  distance  on  the  graph  (~  1  KHZ)  indicating  that  results  of  forcing 
at  frequencies  away  from  natural  frequencies  are  appropriate  for  comparison  to 
the  graph.  With  this  observation  the  results  of  Appendix  A  can  be  combined  to 
indicate  behavior  in  Fig.  4.8.  The  "bending  type"  dissipation  has  the  correct 
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Figure  4.8  Qualitative  malleus  amplitude  vs  frequency  for  two  dissipation 
mechanisms.  Solid  curves  represent  data  from  Manley  and 
Johnstone  (1974)  for  a  guinea  pig  at  nn  R  S^.. 


shape  to  match  the  data,  while  the  "transverse  velocity  type"  dissipation 
indicates  excessive  malleus  motion  at  high  frequencies. 

4.8  Concluding  Remarks 

Some  simple  analogies  have  been  applied  to  address  dominant  mechanisms 
governing  behavior  of  the  tympanic  membrane.  At  100  dB  SPL  membrane  type 
restoring  forces  dominate  over  bending  and  shear  restoring  forces.  At  lower 
sound  pressure  levels,  evidence  indicates  existence  of  a  prestress  maintain¬ 
ing  dominance  of  the  membrane  terms.  Membrane  stresses  clearly  dominate  for 
SPL's  near  or  above  100  dB.  This  observation  indicates  that  a  curvilinear 
membrane  shell  formulation  may  be  sufficient  to  describe  primary  behavior  of 
the  restoring  forces. 

In  addition  to  the  restoring  mechanism,  it  was  determined  that  bending 
type  damping  is  more  reasonable  than  transverse  type  damping.  The  magnitude 
of  damping  may  be  determined  by  comparison  to  malleus  amplitude  vs.  frequency 
curves  or  phase  shift  data. 

Behavior  of  a  string  constrained  by  a  linear  spring  at  one  end  and 
fixed  at  the  other  end  was  used  to  illustrate  concepts  and  ideas.  The  string 
is  not  to  be  misconstrued  as  a  TM  model  -  it  is  presented  simply  as  an  explan 
atory  device.  Readers  interested  in  details  of  the  string  problem  should  see 
Appendix  A. 
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4.9  APPENDIX  4.1:  LINEAR  STRING  PROBLEM 

Consider  vibration  of  a  string  supported  by  a  spring  on  one  end  and 
clamped  on  the  other.  Including  small  bending  and  small  transverse  damping. 


Figure  4.9  Simple  string  subject  to  a  spring  constraint  and  a  fixed  end. 


the  string  motion  is  described  by 


uxx  +  e8  uxxt  +  eTut  ’  \z  Utt  =’TSin  * 


where  the  boundary  conditions  are 

u(0,t)  =  0  , 


and 


ux(*.t)  =  *  u(4,t)  , 


(4.9.1) 


where  T  is  the  tension,  c  is  the  sound  speed,  Pq  is  the  pressure  amplitude, 

a  is  the  forcing  frequency,  and  u  is  the  transverse  displacement.  For  the 

case  of  no  damping  eB  =  ej  =  0,  a  direct  expansion  may  be  applied  to  find  an 

asymptotic  form  for'small  endpoint  motion  ( |_I _ |  <<  1).  Ignoring  transients, 

T-k£ 

the  result  is 


oo 

u(x,t,e)  =  Z  2c2sinut{-Po  Sjn  nitx 
n=l  T(</-</)  nir  l 
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'  i  “  —  ■  — wp  y  • 

n=l  l(w  -to  ) 


Also,  k  =  spring  const.,  Po  =  amplitude  of  applied  pressure,  T  =  aA  =  tension 
c  =  /T/p"  =  sound  speed,  I  =  length,  to  =  forcing  frequency,  =  nirc/£  =  natur 
al  frequency,  1  »  e  =  small  endpoint  motion  parameter. 

Since  the  endpoint  motion  is  small  it  is  advantageous  to  remove  the 
endpoint  motion  from  the  boundary  conditions  and  place  it  in  the  differential 
equation.  To  do  so  let 

(4.9.3) 


u  =  v  -  e  vx(*»t)  x  , 


SO 


v  +  e„v  ^  +  e 
xx  B  xxt 


(vtf  e  =  9(*)  sin  > 


v(0,t)  =  0  , 


where 
and 

V ( £ ,t )  =  0 

For  small  e  we  can  use  a  direct  expansion  for  the  displacement  as  follows: 


(4.9.4) 


v(x,t;e)  =  v0(x,t)  +  e  v^x.t)  +  ...  . 


(4.9.5) 


Using  a  multiscale  analysis  for  <o  =  top  we  find  for  bending  damping  (er  =  0) 

3 

V(B)  a  g*c2  ( i_exp  r-eBto>nt)3  cos  u)  t  sin 
O  T  v  v  2  n  l 

(4.9.6) 


eB  wn 


where 


g*  =  <  -g ( x ) ,  sin  n,|rx  > 

l 


By  a  similar  method  for  to  =  <o  and  eR  =  0,  the  result  for  transverse  damping 


«-*■  **-  ■*- »-  -*-  if  -  ■*  V-  -  -  - >  *-• 


,("0  ~  g*  (i-Pxp("eTc  t]]  cos  a,  t  sin  nirX 
"  v  2  n  £ 


eT“n 


where 


(4.9.7) 


g*  =  <  -g(x) ,  sin  n7rx  >  . 

£ 


If  (1)  *  ID 

n 

the  undamped  case 


then  both  v^T)  and  v^  have  ur2  dependence  very  similar  to 
.  Motion  at  the  spring  is  given  by 


u(£,t)  =  v(£,t)  -  e  v*(i,t)z  , 

and  so 

u(*,t)  =  e  vx(£,t)£  . 


(4.9.8) 


Thus,  the  spring  motion  to  order  e  is  given  directly  as  the  spatial  derivative 


CHAPTER  5 


FORMULATION  OF  A  DYNAMIC  FIBER  COMPOSITE  MEMBRANE  SHELL  MODEL  OF  THE 

TYMPANIC  MEMBRANE 

5.1  Equilbrium  Equations 

The  previous  chapter  uses  the  tympanic  membrane  ultrastructure  and 
experimental  test  data  to  argue  that  membrane  type  stresses  represent  the 
primary  restoring  force  in  the  eardrum.  If  we  neglect  bending  terms  and 
include  only  membrane  stresses,  then  conservation  of  linear  momentum  written 
in  component  form  gives  the  following  curvilinear  membrane  force  equilibrium 
equations  (Novozhilov,  1959). 


1  U_{H  N  )  +  N  +  f(HN  )  -  N  }  +  P  =  0  (5.1.r 

H1H2  3qi  aC|2  3q2  3qi 


1  {  3-(H  N I  )  +  ^  N21  +  ^H2N12^  -  P-  U  +  P2  *  0  (5.1.2! 

H1H2  3q2  3qi  3qi  3q2 


_i_{_i4H2(0iNi  ♦  82ni2))  +  JL(*2\  ♦  6^))} 


(5.1.3: 


.  Nl  .  N2  +  p  0  . 

R  R  3 
I  2 

where  N. .  are  the  products  of  membrane  stresses  with  local  thickness,  R.  are 

I  J  1 

radii  of  curvature,  q^  are  the  curvilinear  coordinates,  P^  are  generalized 
D'Alembert  pressures,  ^  are  coordinate  metrics,  and  ei  are  transverse  defor¬ 
mation  gradients.  As  written,  damping  and  inertia  terms  are  contained  in  P..  ( 

The  force  equilibrium  equations  when  combined  with  an  appropriate  con¬ 
stitutive  law  describe  dynamic  behavior  of  general  curvilinear  membranes.  Fc 
the  eardrum,  coupling  to  the  ossicular  chain  appears  as  a  boundary  condition 
along  the  malleus-tympanic  membrane  attachment.  This  boundary  condition  is 
not  simple  and  in  its  full  form  contains  terms  describing  dynamics  of  the 
ossicular  chain  and  cochlea.  The  boundary  condition  at  the  malleus  couples 


the  membrane  equations  to  the  middle  and  inner  ear.  In  addition  to  the 
malleus  boundary  condition,  the  tympanic  membrane  is  supported  around  its 
outside  periphery  at  the  annular  ring  providing  a  second  boundary  condition. 
Since  we  are  interested  in  describing  both  steady  and  transient  dynamic 
behavior,  initial  displacement  and  velocity  fields  must  also  be  prescribed. 
Any  initial  tympanic  cavity  pressure  or  dynamic  changes  in  the  middle  ear 
pressure  are  contained  in  P3. 

Rather  than  application  of  a  numerical  method,  it  is  desirable  to  find 
an  approximate  closed  form  solution  describing  dynamic  displacement  of  the 
tympanic  membrane.  To  do  so,  we  begin  by  determining  geometrically  dependent 
terms  in  the  equilibrium  equations  for  a  general  eardrum  shape.  The  shape  of 
the  eardrum  can  be  described  as  a  "perturbed  cone"  where  departure  from  the 
perfect  cone  shape  is  given  by  eg(r,0).  The  resting  position  is  defined  by 
f(r,9)  =  ar  +  eg(r,9)  as  shown  in  Figure  5.2.  Definition  of  the  eardrum 


Figure  5.1  Equilibrium  of  an  element  of  composite  material. 
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Figure  5.2  Geometrical  sketch  of  the  tympanic  membrane. 


Selecting  polar  coordinates  qx  =  r  and  q2  *  8  defines  the  position 


vector  r. 


r  *  r  cose  i_  +  r  sine  j,  +  (ar  +  eg)k  , 


where  i,  j,  are  in  the  plane  of  the  annular  ring  and  k  denotes  the  axial 
direction.  On  the  surface,  tangent  vectors  to  coordinate  lines  are 

r*  »  11  cose  i  +  sine  i  +  (a  +  egr)k  (5' 


r2  =  _  =  -  r  sine  i  +  r  cose  j  +  egQ  k  . 

Using  these  the  coordinate  metric  normalizers  are 

H2  =  r  +  0(e2)  , 


H,  *  (l+a2)'/2  *  e°gr  ....  *  0(e2)  , 
2( 1+a2) 


Hi2  =  eage  +  °(e2)  • 


After  finding  unit  normals  to  the  surface,  principal  radii  of  curvature  are 
determined  to  be 


and 


eR  =  ~1  +  e  r(l+a2)  +  o(e2)  , 

1  noQrr  no  a 

R2  =  r(H-a2l  +  o(e)  ,  (5.1.7) 

n0  a 

n  =  (l+a2)-1/2  -  e - +  0(e2)  . 

0  Z(l+a2)s/i 


Also,  letting  w(r,0,t)  denote  the  deformation  normal  to  the  inital  position 
(transverse  deformation),  then 


°i  *  I?  (tl*“2)I/2  - e  «v  *  °<e2>» 


(5.1.8) 


and 

9  =  {I  +  0(e2)}  . 

2  30  r 

The  membrane  stresses  in  the  tympanic  membrane  are  related  to  strains 
and  curvatures  through  a  constitutive  law.  The  most  general  linear  consti¬ 
tutive  law  for  membrane  stresses  in  a  composite  has  the  form  (Jones,  1979) 

N  =  £  e  +  k  ,  (5.1.9) 

where  Ni  is  the  membrane  force,  Q-jj  is  the  membrane  stiffness,  ei  is  the 
strain,  B-jj  represents  coupling  to  changes  in  curvature,  and  ki  is  the  change 
in  curvature.  Since  the  TM  ul trastructure  is  not  symmetric  through  its  cross 
section,  J3  will  be  nonzero  -  a  complication  not  encountered  with  isotropic 
materials.  The  curvatures  are 

k  =  -1  ££l  "  92  3H1  » 

1  TC  3q  H  H  3q 
1^1  1  2  h2 


k 


2 


1  302  -  01  3H2 

TT7  "aqT  hX-sTT  ’ 
2  ^2  12  ^ 1 


(5.1.10 


k  « -  I  (i_l!2  +  . 


and 


Defining  u^r.e)  and  u2(r,0)  as  deformations  tangent  to  the  surface  then 


3  1  e2  +  1  3ul  +  1  3H1  u 


1  2  1  Hx  aqi  H1H2  3q 


(5.1.11) 


and 


e  =  1  e2  +  1  *"2  +  1  3H2  tl 

2  2  2  h2  aq2  h1h2  aqt  1 


'12 


=  ee  +i__i«i .  +  i  aHi  u. 

1  2  3qL  H2  3q2  .  3qi  2  H1H2  3q2 


To  insure  existence  of  a  single  valued  deformation  field  associated  with 
the  stress  forms  we  need  to  satisfy  the  following  compatibility  condition 
(Novozhilov,  1959). 


_J_  [i_  JL(H  ej  -  L.^1  g  -11  _1_( H2y)] 

an  lH  an  2  2  u  an  1  J  7  3d  1 


3qj  Hj  3qi 


Hx  3qt  1  2  Hj  392 


+  _L 


3q, 


[1 _ fL(H  e  )  -  l_i!ll  e  -11  _L.(H2y)] 

H,  3^7  1  H,  3q,  2  2^3q/  *TJJ 


■2  3<2 


=  13  fi  _LJh  02)  -  1  iil  e2  .  l  _i_(H28  0  )] 

T  a7"  Lir  aTT4'  2°2j  rTT'’l  H  2JJ 


2  3q^  Hr  3qx  ^  ^  3qi  1  H‘  3<12 


+  I JL.  [i _ L-(h,92)  -  1  _L.(h20  ej] 

2  3q,  LH,  3q,  11  H2  3^2  ^  9\ 


2  "2  ^2  "2  ^2 
With  this  result,  the  general  formulation  is  complete. 


5.2  Ultrastructure  Simplifications 

The  middle  plane  spacing  of  the  fibers  within  the  tympanic  membrane 
indicates  that  the  bending  -  membrane  coupling  coefficients  B-j j  are  relatively 
small.  For  a  symmetric  laminate  composite  _B  is  identically  zero.  In  addition 
to  small  bending  -  membrane  coupling,  the  membrane  shear  stiffness  Q33  and 
shear  couplings  Q  ,  Q  .  are  small  for  the  eardrum  material.  Neglecting 


these  terms  in  (5.1.9)  gives 
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(5.2.1) 

Moreover,  N12  =  0,  Mi j  *  0,  and  B.  -  0.  If  we  also  neglect  in-plane  inertia 
then  P1  =  0  and  ?2  =  0.  With  these  simplifications  the  equilibrium  equations 
(5.1.1)  -  5.1.3)  reduce  to 

-2  (H  N  )  -Sn  =  0,  (5.2.2) 

3r  2  1  3r  2 

JL  (H  N  )  -  iHl  N ,  =  0  ,  (5.2.3) 

36  k  1  2J  an  1 


36 


and 


L _  {J.  (H  0  N, )  +  (H,  9  N  ) }  -  Hi  -  M  +  P  =  0  .  (5.2.4) 

HxH2  K  2  1  30  ^  1  2  2"  R.  R_  3 


Using  (5.1.6)  the  first  two  equilibrium  equations  can  be  combined  to  show  that 
the  membrane  forces  have  the  form 


N2  =  N2(r)  +  0(e) 

p 

N  =  I  /  N  (n)3n  +  -ElilL  +  0(e) 
1  r  „  2  r 


(5.2.5) 


where  Ci(6)  is  a  constant  of  integration.  Using  this  in  the  third  equilibrium 
equation  (5.2.4)  gives 


J_  {H  (0,  +  _1  ( 0?H. ) )  +  r  -  1 
rH.  N  1  30  2  1  3r  U 


TT 

i 


N2  +  p3  (5.2.6) 

NR  N 
1  2  1 


+  0(e)  =  0  . 

Nondimensionalizing  r  by  the  average  radius  a,  and  w(r,0)  by  the  experimentally 
observed  maximum  transverse  deformation  shows  that  the  ratio  of  the  average 
circumferential  force  "N  to  the  average  radial  force  provides  a  second 
natural  small  parameter  <  =  "N  / N-  .  Representing  w(r,0,t)  by  an  expansion  in 


<  we  have 


5 


w(r,0,t;<)  =  Wo(r,0,t)  +  k  wl(r,d,t)  +  .  .< 

Substitution  into  the  reduced  equilibrium  equations  (5.2.2,  5.2.3,  5.2.6) 
yields  the  following  sequence  of  problems  when  making  use  of  relations  (5.1.8). 
0(1) 

q  ,  S  ( 0)  <  r  <  S  (8)  ,  0  <  e  <  2.  (5.2.8) 

Tr7^  JPTt  ^3 17^  1  2 

where 

wo(S2,0,t)  =  y$[w0] 


wo(S2,9,t)  =  0  . 

y<)>  is  an  operator  describi.a  the  malleus  boundary  condition  including  ossicular 
and  cochlear  dynamic  terms.  Assuming  e/<  =  0(1),  then  the  order  <  problem  on 
the  same  domain  is 
0(<) 

32wi  +  ^  33wi  =  i  32wi  r  i  32w0  .  c  33w0 

TF~  3r*3t  t7  Tt7~  r  3r3t  "FlFTt 

,  (5.2.9) 

-  {_!_ +  (l+<*2)  a  wo  -  (l+a2)rl 
rNL  3r  r2Nt  302  R2 

+  £  {I_  +  £_  +  JL  a2wo} 

<2  Rx  N1  It1' 

where 

^(Sj.e.t)  =  o  ,  O 


and 


w.(S,,0,t)  =  0  , 


where  a  D'Alembert  reverse  inertial  pressure  has  been  included  transverse 
to  the  resting  position.  The  boundary  r  =  S^e)  represents  projection  of 
the  malleus  boundary  onto  the  (r,9)  plane,  and  r  =  S2(9)  represents  pro¬ 
jection  of  the  annular  ring  boundary  onto  the  (r,9)  plane.  The  c  terms 
represent  first  order  linear  bending  type  damping  (see  section  4.1). 

If  the  radial  membrane  force  N1  were  known  then  the  0(1)  problem 
could  be  solved  directly. 

As  e  *  0  and  the  shell  becomes  a  cone  we  find  N1  ♦  aN1(9)/r.  It  can 
also  be  shown  by  looking  at  equilibrium  equation  (5.2.2)  that  N1  *  aN^eJ/r 
as  k  ♦  0.  Hence  it  seems  very  reasonable  to  let 


N  =  aNl<9) 
1  r 


(5.2.10 


for  small  k  and  small  e.  Using  this  form  it  is  straightforward  to  estimate 
N1  using  the  orthotropic  constitutive  law  approximations  of  Chapter  3. 

The  form  of  equation  (5.2.8)  suggests  that  it  may  be  appropriate  to 
model  the  tympanic  membrane  as  a  set  of  coupled  radial  fibers.  This 
alternative  approach  is  outlined  in  Appendix  A. 


5.3  Asymptotic  Solution  For  w0 

Consider  the  case  when  the  damping  and  sound  speed  are  constant  to 
first  order  such  that  c  and  c  appear  as  constant  coefficients.  To  homogenize 
the  boundary  conditions  in  (5.2.8)  introduce  a  transformation  on  the  depen¬ 


dent  variable 


W„(r,e,t)  -  u  (r,9,t)  *  y  ♦[»„] 


(5.3.1) 


The  first  order  problem  becomes 


^  *  5  -!l2>  -  ‘  (1^°  *  1^2  Y  i!±)  ♦  q 

ir*  ir^n  C7  9tz  s,-s,  at7 


(5.3.2) 


-  75  - 


with  homogeneous  boundary  conditions 


Uq  (Sx ,e,t)  =  0  , 


u0  (S2,e,t)  =  0 


where  y  is  a  third  small  parameter  representing  small  motion  of  the  malleus 
and  <(>  is  an  operator  describing  middle  ear  and  cochlear  connections. 

Looking  for  an  asymptotic  solution  for  small  y  ,  consider  a  direct 
expansion  of  uo(r,0,t) 


u0(r,9,t;  )  =  v0(r , 0 ,t)  +  y  v^r.e.t)  + 


(5.3.3) 


Substitution  of  this  into  (5.3.2)  leads  to  the  following  sequence  of  problems 


32v0  +  c  33 v0 

JF~  ar2at 


where 


+  SL(0)  <  r  <  $2(0)  ,  0  <  0  <  2n  (5.3.4) 


v0(S1,9,t)  =  0 


vo(S2,0,t)  =  0 


a2vi  +  c  a3vi  =  I  a2vi  +  r-S2  324>[  V0  ] 

3 r2  3r23t  C^Tt7""  C*(S,-S,)  Tt7 


(5.3.5) 


v^S^O.t)  =  0 


vl(S2,0,t)  =  0  . 


where 


For  the  case  of  harmonic  forcing  q  =  qo(r,0)  sin  wt,  the  solution  to  (5.3.4) 
is  easily  found  by  eigenfunction  expansion  to  be 


»  c2  qn  sin(uit  +  i »n) 
vo  ■  £  *' 


n=l  ulJ[  l-iO-W  +  U*)2- 

n  V  ain 


(5.3.6) 
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where  c  is  the  damping  coefficient,  c  is  the  sound  speed,  cn,  ^  are  constants 
determined  from  the  inital  conditions, 

/  2  \ 

if/  (0)  =  Arctan  f^n 
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(5.3.7) 


q  (9) 
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2_-  /  1  q0(r , e)  sin  )J  dr. 
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^2 


To  solve  the  0(y)  problem  in  (5.3.5)  it  is  necessary  to  find  the  malleus 
deformation.  To  determine  the  malleus  motion>  define  the  differential  operator 
L  relating  angular  deformation  of  the  malleus  <t>  to  the  moment  transmitted 
through  the  tympanic  membrane  (see  Fig.  5.2). 


L  [♦]  =  M(t) 


(5.3.8) 


Expanding  the  moment  M  =  Mq  +  y  +  0(y2)  allows  direct  calculation  of  the 
change  in  moment  A  as  an  integral  of  forces  on  the  malleus  imparted  by  tym¬ 
panic  membrane  fibers. 


A  M  (t)  =  A  /  (N .  x  0) *m  ds 
u  c  1 


(5.3.9) 


Where  is  the  tension  vector  in  the  direction  of  a  radial  fiber  at  the 

malleus,  0  is  the  position  vector  from  the  malleus  instant  center  of  velocity 


to  Sj ,  and  m  is  a  unit  vector  along  the  instantaneous  axis  of  rotation  of  the 

malleus  (see  Appendix  B).  This  gives  M  explicitly  in  terms  of  V  allowing 

o  o 

direct  calculation  of  the  malleus  motion  by  solving  the  differential  system 

A 

♦(t)  =  L“1[M0(t) J  (5.3.10) 

From  the  malleus  angular  deformation  we  can  determine  the  deformation  at  a 
particular  point  on  the  malleus 

^  A  A  A 

y  <t>0  [v0  ]  s  (0  x  ♦  (5.3.11) 

A 

where  $  is  the  instantaneous  rotation  of  the  malleus  and  s  is  the  unit  normal 
in  the  direction  of  the  component  being  calculated.  With  this  result,  the 
0(y)  problem  has  the  same  form  as  in  (5.3.6)  where  q  is  replaced  by  q*.  given 
as 

(5.3.12) 

Knowing  <t>Q ,  VQ  and  determines  wQ  to  0(y2). 

w0(r,e,t)  =  VQ(r,e,t)  +  y  V  (r.e.t)  +  y  *q  +  0(y2)  .  (5.3.13) 

Sl'S2 

The  above  expression  is  a  closed  form  asymptotic  solution  to  the  general  tym¬ 
panic  membrane  problem.  It  is  coupled  to  the  ossicles  by  membrane  stresses 
and  includes  effects  of  middle  ear  and  cochlear  dynamics  through  the  very 
general  operator  L.  The  outer  ear  and  middle  ear  air  chambers  are  indirectly 
included  as  perturbations  on  the  excitation  pressure  field. 

5.5  Qualitative  Numerical  Results 

Representing  the  geometry  of  a  cat  we  selected  an  annular  ring  superior 


radius  of  0.35  cm,  and  anterior  radius  of  0.30  cm,  malleus  angle  defining 
depth  of  the  drum  of  0.64  radians,  average  density  of  2.0  grams/cm3, 
average  thickness  of  40  x  10-4  cm,  and  nondimensional  damping  coefficient  of 
3  x  10-5.  In  addition  to  these  parameters  the  local  stiffnesses  were  esti¬ 
mated  using  equation  (3.6.11)  along  with  ul trastructural  observations  of  the 
volume  fractions  of  radial  and  circumferential  fibers.  A  very  small  initial 
prestress  was  imposed  by  rotating  the  malleus  10-3  radians  with  respect  to  its 
no  force  position.  This  initial  rotation  could  easily  be  imparted  by  the 
tensor  tympani  or  ligament  preloads. 

Using  these  values,  equation  (5.3.13)  reproduces  vibrational  shapes 
observed  experimentally.  Figure  5.3  indicates  vibration  of  the  drum  when 
subjected  to  a  pure  600  Hz  tone  at  105  dB  SPL.  The  vibrational  shape  essen¬ 
tially  matches  that  observed  experimentally  by  Khanna  and  Tonndorf  (1972). 

The  shapes  shown  in  Figure  5.3  are  based  on  a  simplified  geometry  and  neglect 
the  pars  flaccida  and  hence  the  shapes  deviate  slightly  from  the  results  of 
Khanna  and  Tonndorf.  The  asymptotic  solution  given  by  (5.3.13)  is  capable  of 
describing  behavior  for  the  more  complex  geometry  and  material,  however  the 
intent  here  is  simply  to  look  at  qualitative  behavior  leaving  detailed 
calculations  for  Part  II  of  this  report. 

At  600  Hz  the  drum  vibrates  in  its  fundamental  shape  essentially  in 
phase  with  the  excitation  pressure.  As  the  frequency  is  increased  the  mem¬ 
brane  vibrational  amplitude  decreases  and  begins  to  lag  the  excitation.  Above 
4000  Hz  a  transition  from  the  simple  fundamental  vibrational  shape' to  a 
complex  shape  is  reproduced  by  the  asymptotic  solution.  When  forced  at  higher 
frequencies,  the  vibrational  shape  of  the  tympanic  membrane  is  constantly 
changing  in  time.  Figure  5.4  shows  normalized  displacement  at  4000  Hz  for 
two  different  times  during  the  excitation  cycle.  Time  averaged  experimental 


600  Hz 


6000  Hz 


results  represent  a  weighted  average  of  many  curves,  as  shown  in  Fig.  5.4, 
over  numerous  cycles  (time  averaged  holography  for  example). 

Results  that  can  be  obtained  from  the  asymptotic  solution  are  numerous 
The  effects  of  geometrical  changes,  material  alterations,  scar  tissue,  and 
middle  ear  changes  represent  only  a  few  of  the  possible  topics.  In  addition 
to  these,  many  clinical  applications  may  be  found.  Spatial  integrations  of 
the  asymptotic  solution  can  be  compared  directly  to  tympanograms  rnd  acoustic 
impedance  data.  The  comparisons  can  be  used  to  make  the  tympanic  membrane 
transparent  to  clinicians  or  to  assess  function  of  the  eardrum  itself.  Part 
II  of  this  report  will  study  the  asymptotic  solution  and  address  some  of  the 
topics  mentioned. 


5.5  APPENDIX  A 
A  FINITE  RADIAL  FIBER  MODEL 

By  studying  the  behavior  of  the  TM  under  various  conditions,  it 
appears  that  the  radial  (or  meridian)  fibers  play  the  dominant  role  in  the 
function  of  the  drum.  Any  axial  stresses  in  the  meridian  fibers  act 
directly  on  the  malleus,  near  the  umbo  in  most  cases,  and  thus  transmit 
force  directly  to  the  ossicles.  The  complicated  vibrational  shape  at  high 
frequencies  observed  experimentally  also  suggests  that  the  circumferential 
fibers  play  a  secondary  coupling  role  while  the  radial  fibers  dominate  the 
shape  and  transmittance  behavior.  Assuming  that  this  postulated  radial 
dominance  exists,  it  is  possible  to  formulate  a  coupled  set  of  m  partial 
differential  equations,  that  individually  represent  the  physics  of  a  sec¬ 
tional  strip  of  the  TM.  Since  most  of  the  radial  fibers  run  in  the  meridian 
direction  originating  at  the  umbo,  most  of  the  strips  will  be  pie  shaped 
sections  of  the  membrane.  Each  radial  structure,  or  pie  shaped  strip,  is 
attached  to  the  annular  ring  on  the  outside  edge  and  attached  to  the  malleus 
on  the  inside  edge.  It  is  assumed  that  the  fiber  elements  interact  with 
each  other  through  a  membrane  type  tension  term,  however  the  element  itself 
contains  both  bending  and  tension  energy  storage  mechanisms.  Dissipation  of 
energy  is  assumed  to  be  of  the  bending  type  and  hence  will  be  proportional 
to  time  rate  of  changes  in  the  two  principle  radii  of  curvature  of  the  drum 
surface. 

We  begin  the  model  by  defining  N  coordinates  X^  originating  at  the 
malleus  attachment  of  each  fiber  directed  outward  toward  the  fiber  attach¬ 
ment  at  the  annular  ring.  The  X^  coordinates  are  along  straight  lines,  such 
that  the  middle  plane  of  each  radial  element  in  its  original  static 
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equilibrium  position  is  defined  by  Uqn(Xn) .  For  fibers  attached  at  the 
umbo,  the  coordinate  will  correspond  to  the  polar  coordinate  r;  for 
fibers  attached  on  the  shank  of  the  malleus,  the  coordinate  will  cut 
across  polar  coordinate  lines. 

In  order  to  model  one  particular  element,  imagine  that  we  have  cut 
the  fiber  and  its  surrounding  tissue  from  the  membrane.  The  single  radial 
fiber  (or  group  of  radial  fibers)  imbedded  in  the  tissue  slice  is  under  a 
meridian  stress  amn,  has  thickness  hN(Xf|),  and  width  wn(Xn).  The  static 
equilibrium  position  of  the  Nth  fiber  element  is  described  in  global 
coordinates  as  z*  =  f(r,8),  or  in  local  element  coordinates  as  Uon(Xn)* 
Figure  5.5  indicates  the  coordinates  for  a  typical  radial  element. 

Viewing  the  radial  element  from  the  side,  it  appears  as  a  curved 
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Figure  5.5  Sketch  of  a  pie  shaped  slice  of  the  tympanic  membrane. 
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beam  under  axial  tension  forced  by  the  applied  pressure  and  adjacent  ele¬ 
ments.  The  fiber  element,  if  originating  from  the  umbo,  will  have  a 
variable  width  nearly  proportional  to  the  lengthwise  coordinate  (or  r). 
Denote  the  deflection  of  a  point  perpendicular  to  the  static  equilibrium 
position  as  UN(XN,t),  and  the  component  of  deflection  normal  to  UN(XN,t) 
as  V|\|(XN,t).  Displacement  at  the  malleus  attachment  is  denoted  (UNE(0,t), 
VNE(0,t)). 


Figure  5.6  Side  view  of  a  single  radial  fiber  in  local  coordinates. 

Suppose  the  element  thickness  w  <<  a  (a=length)  and  that  the  initial 
curvature  is  much  larger  than  the  thickness  (or  in  fact  the  length),  such 
that  a  linear  stress  distribution  due  to  internal  bending  moment(sl  is  a 
reasonable  approximation.  For  the  same  reason  it  is  assumed  that  the  axial 
stress  distribution  can  be  approximated  as  uniform  across  each  section. 


Under  such  conditions  the  primary  forces  acting  on  a  differential  slice  of 
the  radial  element  in  X  are  indicated  in  Fig.  5.7. 


Figure  5.7  Equilibrium  of  an  element  of  drum  material. 

Considering  conservation  of  momentum  normal  to  the  membrane  surface 
we  can  derive  a  differential  equation  describing  the  element.  Noting  that 
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Expanding  for  small  deformation,  this  can  be  written  as 
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The  previous  equation  reduces  to  the  first  order  membrane  shell  model 
as  the  bending  stiffness  approaches  zero.  Equation  (5.5.3)  is  not  exact,  but 
does  contain  bending  terms  and  circumferential  coupling  not  contained  in  the 
first  order  shell  membrane  equation. 
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APPENDIX  B 

MALLEUS  MOMENT  AND  DISPLACEMENT  RELATIONSHIPS 

In  order  to  determine  the  moment  imposed  by  the  tympanic  membrane  on 
the  malleus  about  its  instant  center  of  velocity,  the  radial  tension  vector 
T  and  position  0  are  defined  as  indicated  in  the  diagram  below 


i: 

t. 


Figure  5.8  Sketch  indicating  tension  vector  due  to  a  radial  fiber  on  the 
malleus 

The  tension  vector  is 
♦  . 

T  =  Tcos  a  (cos8  i  +  sine  j]  +  Tsin  a  k 

J  A  A  A  A 

=  - [cose  i  +  sine  a  +  a  k[  (5.5.^ 


=  t  i  +  t  j  +  t  k 
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and  the  position  vector  is 
+• 

A  A  A 

0  =  (i-Bcose)(-i )  +  Bsine  j  -  (t-Bcose)  sin  k 
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The  resulting  moment  is 
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The  component  causing  rotation  about  the  instant  center  of  velocity  is 
M  =  M  .  m 
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where  we  have  taken  the  instant  centerline  to  be  in  the  plane  of  the  annular 

ring.  This  assumption  could  easily  be  relaxed  by  including  a  third  component 
in  m. 

During  a  deformation  a  and  T  change.  To  obtain  the  change  in  the 
moment  during  the  deformation  let 


a(t)  =  a  +  Aa 
T(t)  =  T  +  AT 


(5.5.8) 


Substituting  into  the  moment  expression  and  using  a  »  Aa  ,  T  »  AT  we  obtain 
the  change  in  moment 
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The  displacement  of  the  malleus  at  the  attachment  of  a  radial  fiber 
is  defined  by  vector  3  and  is  related  to  malleus  rotation  <j>  by 
+  +  *  - 

d  =  (0  x  M) <J>  (5.5.10) 

^  A  A  A 
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The  component  in  the  k  direction  is 

d  =  {(ii-Bcose)  cos  8  -  B$inesin6}4> 
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APPENDIX  C 
IMPULSE  RESPONSE 

If  the  tympanic  membrane  is  subjected  to  an  impulse  of  magnitude  I  at 
time  t  =  0  then  the  V0  function  is  defined  by 


92V0  .  t  33V0  .  1  32»0  ,  :  J(t) 
3r2  3r23t  C7  It7" 


(5.5.12) 


For  an  eardrum  initially  at  rest  Vn(r,0,O)  =  0  and  Vg  (r,0,O)  =  0.  Integrating 


(5.5.12)  across  the  impulse  in  time  allows  us  to  reformulate  the  problem  as  a 
homogeneous  operator  with  nonzero  initial  conditions.  The  equivalent  problem 
i  s 


32V0  +  33Vq  _  i  32V0  ,  o 


(5.5.13) 
.2, 


Sr2  "  ar^at  V-  3t2 
with  initial  conditions  at  time  zero  plus  VQ(r,0,O+)  =  0  and  VQ(r,0,  0+)=  -Cgl. 
The  general  solution  to  (5.5.13)  can  be  written 
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Imposing  boundary  conditions  we  find 
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C  may  be  found  by  taking  the  inner  product  of  (5.5.15)  with  sin  (nir  Sl—l . ) 
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For  constant  I  we  find 
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CHAPTER  6 

A  CONSISTENT  MATHEMATICAL  MODEL  OF  THE  EUSTACHIAN  TUBE  AND 
MIDDLE-EAR  AIR  CHAMBERS 

6.1  Importance  of  Wave  Motion  in  the  Middle  Ear  Air  Chamber(s) 

In  order  to  formulate  a  simple  model  of  the  middle  ear  air  chamber(s) 
that  adequately  describes  action  of  the  chambers,  we  must  determine  if  it  is 
reasonable  to  neglect  wave  motion  within  the  chambers.  As  a  rough  estimate  of 
the  importance  of  wave  motion,  we  may  compare  one  quarter  wavelength  of  the 
forcing  pressure  to  a  length  characterizing  the  primary  air  chamber.  Assuming 
that  the  sound  speeds  inside  the  chambers  and  outside  the  external  ear  have 
similar  values,  then  the  wavelength  comparison  amounts  to  comparing  the  period 
of  the  forcing  signal  to  the  time  required  for  a  signal  to  traverse  the  middle 
ear  air  chamber.  If  the  traversing  time  is  much  less  than  the  forcing  period, 
then  fluid  properties  within  the  primary  air  chamber  will  be  essentially 
uniform  and  wave  propagation  can  be  ignored.  On  the  other  hand,  if  the 
traversing  time  is  of  the  same  order  or  greater  than  the  forcing  period,  then 
wave  propagation  within  the  middle  ear  air  chamber  must  be  considered. 

As  a  conservative  estimate,  the  time  required  for  a  signal  to  traverse 
the  primary  chamber  can  be  estimated  as  the  distance  divided  by  the  sound 
speed.  If  the  characteristic  length  of  the  primary  chamber  is  taken  to  be 

1.6  cm,  and  the  sound  speed  is  approximated  as  350  m/s,  the  traverse  time  is 

4.6  x  io-5  sec  corresponding  to  a  full  wave  frequency  of  22  KHz  and  a  quarter 
wave  resonance  of  5.5  KHz.  Hence,  if  wave  motion  in  the  middle  ear  air 
chamber(s)  is  neglected,  then  the  forcing  frequencies  must  be  limited  to 
values  well  below  5.5  KHz  for  the  example  chamber  size  considered  here. 


6.2  Discrete  Middle  Ear  Air  Chamber  and  Eustachian  Tube  Model 

Suppose  we  model  the  middle  ear  air  cavities  as  consisting  of  one 
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primary  cavity,  connected  to  the  Eustachian  tube,  and  “n"  secondary  cavities 
as  indicated  below. 


Figure  6.1.  Schematic  diaoram  of  air  chambers  and  Eustachian  tube. 

Displacement  of  the  tympanic  membrane  is  denoted  by  w^.  Fluid  velocity 
down  the  Eustachian  tube  at  the  connection  to  the  primary  cavity  is  denoted  by 
Ue,  while  fluid  velocity  into  the  secondary  cavities  is  denoted  Ui  (1*2, . . .  ,n) . 
The  volume  of  each  chamber  is  considered  fixed  for  development  of  this  discrete 
model.  The  fluid  state  is  defined  by  pressure,  Pn  and  density,  pn  in  each 
chamber  volume  Vn. 

As  a  simple  model,  consider  the  case  of  "slow"  TM  vibrations  such  that 
wave  motion  within  the  middle  ear  chambers  can  be  neglected  arrd  the  state  within 
each  chamber  is  considered  uniform, 
a)  Primary  chamber 

Conservation  of  mass  may  be  written  for  a  control  volume  enclosing  the 
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primary  chamber,  traversing  the  connecting  openings  at  the  periphery  of  the 
primary  cavity,  as  follows:  (Potter  and  Foss,  1975) 

f  i£.dV,+  /  p,  l!!I  dA  +  /  p,u  dA.  (6.2.1) 

V2  at  1  TM  1  at  e  1  e  1 


n 

+  Z  /  p  u.  dA  *  0  , 
i=2  A  1 

i 

where  A^  is  the  cross  section  area  of  the  connecting  tubes  (Fig.  6.1).  For 
slow  vibrations  the  properties  within  each  cavity  are  considered  uniform,  and 

the  velocities  at  the  openings  are  approximated  by  their  average  values  such 
that  (6.2.1)  may  be  written 

dpi  +  pi  (/  jjvq-  dA  +  u  A  +  Z  u.  A. )  =  0  .  (6.2.2 

dt  Vx  at  e  e  i,2  1  1 

Suppose  changes  in  density  are  small  and  denoted  by  6(t) ,  then  conservation 
of  mass  for  the  primary  chamber  becomes 

V  po  +  6i(t)  • 

Combining  this  with  (6.2.2)  we  find 


d6x  +  Pofj  jndA  +  uA  +  J  u.  A. )  *  0  .  (6.2.3 

W~  Vj  Hm  at  e  e  i-z  1  1 

b)  Secondary  cavities 

Conservation  of  mass  for  one  of  the  secondary  chambers  may  be  written 

J  dpdV  -  /  p  u  dA  =  0  .  (6.2.4 

V  dt  i  A  n  n 


Again  assuming  uniform  properties  this  reduces  to 

dPn  -  Pn  u 
dt  Vn  n 
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Applying  conservation  of  linear  momentum  to  a  slug  of  fluid  connecting 


the  primary  chamber  to  a  secondary  chamber  we  have 

dun  _  cn  u  =  fp  _  p  1  An  n  =  2,3  ... 
dt  mp  n  1  n  mn 


(6.2.6) 


where  M  =  mass  =  p  A  £  ,  C  =  viscous  damping  and  A  =  cross  sectional  area, 
n  n  n  n  n 

Combine  this  with  (6.2.4)  to  obtain 


3t  Pn  dt  A  £  p  dt  1  n  Pn*n 
n  n  n 


(6.2.7) 


Assuming  isentropic  fluid  behavior  then 

p  =  P-°-  =  constant 
pY  pj 


(6.2.8) 


Also  assuming  changes  in  density  are  small  relative  to  initial  density  p 


define 


pn(t)  a  po  +  V*)  *  for  6n  «  po 


(6.2.9) 


Combining  results,  the  first  order  problem  for  the  perturbations  in  the  density 
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Hence  the  density  in  the  Nth  chamber  may  be  described  by 


■6n  .  C 


Note  that  the  density  is  coupled  to  the  primary  chamber  through  6 


6.3  Eustachian  Tube 


Applying  conservation  of  momentum  to  the  Eustachian  tube  in  the  same 
way  that  it  was  done  for  the  secondary  air  cavities  we  find 

Ji4“es  (pr  PJ  —  * 

dt  mf  1  e 

There  are  two  cases  of  interest;  when  the  tube  is  open  and  when  the  tube  is 
closed. 


Pl'  pl*  chamber  state 


w 


Figure  6.2  Schematic  of  the  Eustachian  tube. 


i )  Open  tube. 

If  the  tube  is  open  then 

V  %  ~  •’)  • 

0 


Combining  this  with  (6.3.1)  gives' 

^ue  +  £e  ue  *  po  Y  t  open  (6.3, 

<Tt”  n»f  p0  mf 

ii)  Closed  tube. 

Apply  conservation  of  mass  to  the  closed  tube  to  obtain  another 


secondary  chamber  equation  in  the  following  form 

d2ge  +  Ce  =  Pq  t  fs  -  6  ) 

dt2  mf  dt  p^-1mf  1  e 

o 

(6.3.3) 

u  A  =  l!i  ,  CLOSED 
e  e  pe  dt 


6.4  Primary  Middle  Ear  Air  Chamber  Acoustic  Model 

If  the  forcing  frequency  is  of  the  order  5  KHz  or  greater,  then  wave 
propagation  within  the  primary  middle  ear  air  chamber  should  be  considered. 
Since  the  secondary  chambers  are  usually  much  smaller  than  the  primary  chamber, 
the  discrete  model  of  their  behavior  will  generally  remain  reasonable  through¬ 
out  the  audible  frequency  range.  Hence,  in  order  to  extend  the  middle  ear  air 
chamber  model  to  be  valid  at  frequencies  above  5,000  Hz  we  must  derive  a  new 
primary  chamber  equation  that  allows  variable  fluid  properties  within  the 
chamber. 

Equation  (6.2.1)  for  the  primary  cavity  remains  valid  and  may  be  written 
as 

a  n 

/  !£.dV+p  Jf  Jzl  dA  +  u  A  +  S  A.  u. )  =  0  .  (6.4.1) 

VL  at  1  °'tm  at  e  6  i=2  1  1 

Introducing  the  velocity  potential  $  within  the  primary  cavity,  acoustic 
perturbations  will  satisfy  the  wave  equation 


The  v  of  the  air  is 


1 


(6.4.2) 


(6.4.3) 


Assuming  adiabadic  motion. 


(6.4.5) 


p  =  -p0(J£)  po  Y 


Substituting  (6.4.5)  into  (6.4.2)  we  have 


1 -y  1-y 

.§£.  »  -p  1  C2?2  <*(.££)  Y  p  Y  > 

at  p°y  o  o 


and  since 


1-Y  1-Y 


p  Y  (2i)  Y  ~  J_  , 

0  M  P0 


then 


i£  ~  JL°  C?72*  • 
at  yP^  0 


Using  Green's  formula  we  find 


2 

Co  -2, 


/  i2~Vz«.dV  +  /  i!!idA  +  uA  +  E  A.  u. 

Vj  TP0  1  TM  at  e  e  i.2  1  1 


and 


/  {£l  n  +  v  }dS  =  0  • 

cs  yp„  s 


(6.4.6) 


(6.4.7) 


(6.4.8) 


s  0  *  (6.4.9) 


(6.4.9) 


But  this  expression  must  be  valid  for  arbitrary  control  surface  selections,  so 
the  boundary  condition  on  the  control  surface  is  just 


1  7$,n  +  v  =0 

p»  s 


(6.4.10] 


Surface 

The  equations  presented  assume  that  small  perturbation  assumptions  of 
acoustics  are  valid  and  complete  the  model. 


6.5  Comments  on  the  Lightly  Forced  or  Damped  Linear  Problem 

Suppose  behavior  of  the  system  can  be  described  by  linear  operators  in 
time  and  space  such  that  the  following  general  form  is  obtained; 

L[u(x)]  +  XM[u(x)]  =  e  F[x,u(n)] 


(6.5.1 


where  x  Is  the  independent  variable,  u  is  the  dependent  variable,  e  is  a  small 
parameter,  and  L,  M  and  F  are  linear  operators.  It  is  assumed  that  the 
boundary  conditions  are  natural,  and  the  general  problem  is  self  adjoint.  To 
avoid  an  order  singularity  we  also  require  n  to  be  less  than  the  highest  order 
derivative  appearing  in  operators  L  or  H  ((n)  denotes  differentiation  n  times). 
The  above  form  is  consistent  with  separated  linear  plate,  shell,  or  membrane 
equations  subject  to  small  damping  or  forcing,  and  hence  is  expected  to  be 
reasonable  for  the  tympanic  membrane  problem. 

Suppose  we  consider  an  asymptotic  expansion  of  the  form 

u(x;e)  *  uQ(x)  +  eu^x)  +  e2u2(x)  +  ...  »  (6.5.2) 


and  let 


A  *  +  eX1  +  e2A2  +  ...  .  (6.5.3) 

Substitution  into  the  original  equation  gives  the  following  sequence  of  problems. 

0(1) 

L[u]  +  Ao  M[u]  -  0  .  (6.5.4) 

Solution  of  the  first  order  problem  will  yield  a  set  of  eigenvalues 
A€[x0n]  and  associated  eigenfunctions  [ ipon ] •  Suppose  there  exists  an 
eigenfunction  expansion  theorem  for  this  problem  such  that  the  general  solution 
may  be  written  as  an  infinite  sum  of  the  eigenfunctions. 


u  (*)  =  E  An'Mx*  *on)  (6.5.5) 

0  n=l 

If  the  original  operator  were  obtained  using  some  sort  of  separation  of  vari¬ 
ables  then  the  "A"  coefficients  will  be  general  functions  dependent  on  the  other 
independent  variables  in  the  problem.  The  exact  form  of  the  functions  can  only 
be  determined  by  Inspecting  the  problem  containing  all  independent  variables. 

An  example  problem  to  follow  will  Illustrate  this. 


The  eigenfunction  expansion  above  will  automatically  meet  all  boundary  condi¬ 
tions  of  the  problem.  Associated  with  the  expansion  it  is  assumed  that  there 
exists  an  inner  product,  such  that  evaluation  of  the  coefficients  is  direct. 
The  next  order  problem  can  be  written 

0(c) 

lcux3  +  yi[ui]  =  +  F£x»uSn)]  •  (6*5,6) 

As  before  using  the  eigenfunction  expansion  theorem,  completeness  of  the  set 
On]  allows  us  to  write  the  solution  of  (6.5.6)  as  a  sum  of  the  eigenfunctions 
as  follows 

=  jj  Wx:  Xon>  •  <6-5-7> 

Substituting  into  the  above  differential  equation 

oo 

E  {BnaOn]  +  XqMC*  ])}  (6.5.8) 

n=l 

=  Z  AnMC4»n3  +  F[x,uJ")]  .  (6-5-9) 

n=l 

For  a  particular  mode  of  the  order  one  problem  we  know 

jj  Bn  *  Xon  *  0 

Hence  we  must  have 

Xln  Z  AnM[<(»n]  =  F[x,ujn)]  =  F[uJn)3  (6.5.9) 

n=l 

Assuming  F  is  linear,  it  can  be  expanded  in  terms  of  the  eigenfunctions  as 

00 

F[u(n)]  =  z  A  FO  ]  (6.5.10] 

0  n=l  n  n 

Combining  the  two  expressions  gives 


X,„  I  AM[t]  »  Z  A_F[*  ] 


(6.5.11 


Taking  the  inner  product  of  each  side  (using  appropriate  weight  for  the 


operator) 


z  An  PC^n]> 

ln  An,  ML4»nJ> 


(6.5.12) 


If  the  problem  is  self  adjoint  then  the  eigenvalue  perturbation  may  be 


written 


x  _  <^»  FC%]> 
lm  <^* 


(6.5.13) 


Using  this,  eigenvalues  of  the  full  problem  to  0(e)  are 


\  =  X  +  e  F,^j>  +  0(e2) 
m  om  <4*m‘  Mt%]> 


(6.5.14) 


Simple  Example  Problem 

Consider  free  vibration  of  a  circular,  isotropic,  flat  membrane  that  is 
fitted  over  a  closed  chamber  filled  with  gas. 


V  =  closed  volume 


Figure  6.3  Circular  membrane  over  a  closed  air  space. 


Assume  that  we  are  interested  in  low  frequency  vibrations  of  the  membrane  such 
that  wave  motion  in  the  air  space  can  be  neglected.  As  explained  earlier  this 
assumption  is  valiji  if  the  time  required  for  a  wave  to  traverse  the  air  space 
is  much  less  than  the  period  of  vibration  of  the  membrane. 

If  damping  is  ignored  then  vibration  of  the  membrane  may  be  approxi¬ 


mated  by 


V2u  -  B  //  udA  =  i_  u.  ,  0  <  r  <  a  ,  0  <  9  <  2ir 


(6.5.15) 


C2pVq  D 


u(a,9,t)  =  0  and  u(O,0,t)  is  bounded  . 

Parameters  in  (6.5.15)  are  B  =  bulk  modulus,  p  =  membrane  density,  VQ  =  initial 
gas  volume,  and  C  =  membrane  sound  speed.  Separating  variables  the  radial 
function  is  of  the  form 


where 


u(r,e,t)  =  A„,  sin  (me  +  <*,)  sin  (u>Kt  +  tc,c)w(r)  , 


<j2w  +  idw+(02-Ei)w  =  B  t  W2irrdr 
dr7  r  dr  r7  T ~  J 

0  0 


(6.5.15) 


If  e  «  1  then  the  previous  analysis  applies.  Writing  the  radial  function  as 


a  straightforward  asymptotic  series 


w  =  wQ  +  ewt  +  e2w2  +  ... 


(6.5.16) 


d2Wo  +  1  dw0  +  (-  2  .  m2]w  = 


dr2  r  dr 


0  , 


wQ(0)  bounded  , 


wQ(a)  =  0  . 


(6.5.17) 


The  solution  to  this  problem  may  be  written  in  terms  of  Bessel  functions 


w .  =  C*  J  (6.  r)  +  C**  Y  (s  r) 

Onm  nm  m  '  Qnm  '  rim  m  '■onm  ‘ 


(6.5.18) 


where  the  second  set  of  constants  must  vanish  to  meet  the  deflection 

condition  at  the  origin,  and  the  eigenvalues  6  are  determined  by  imposing 

onm 

the  second  boundary  condition  at  the  outside  ring. 


The  next  order  problem  is 
0(e) 

L[w1 3  +  82  M[Wl]  =  -82  m[wq]  +  F[w03 

I  _  r2  d2  +  r  d  _  m2  (6.5.20) 

IF*  IF 

M  =  r2 

a 

F  =  r2  /  w„2irrdr 
J  o 
o 

Noting  that  the  left  side  is  zero  for  some  eigenvalue,  perturbation  of  the 
eigenvalues  is  given  by 


/  W0nnmr)r2  /  {Jje0nmn)2nndn}dr 


6inm 


''m'‘"0nm'  r  1  1  nr  Onm 
o  o 


(6.5.21) 


I  r2(Jm(60nn,r))2'lr 
0 


The  eigenvalues  of  the  full  problem  may  be  written  to  order  epsilon  as 


B2  »  3? _  +  e  B2  +  Q( e2) 


(6.5.22) 
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